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ABSTRACT 



ANGULAR DISTRIBUTION OF J/PSI DECAYS 
IN DIMUON CHANNEL IN 800 GEV 
PROTON-COPPER COLLISIONS 
BY 

TING-HUA CHANG 

Doctor of Philosophy in Physics 
New Mexico State University 
Las Cruces, New Mexico, 1999 
Dr. Vassih Papavassihou, Chair 

The angular distribution of decays in the channel in 800 GeV 

proton-copper coUisions has been measured for xp > 0.25. The polarization pa- 
rameter A is extracted in 1 GeV of pt and 0.1 oi xp bins for two magnet configu- 
rations with different acceptances. The data indicate that the J/ip^s are produced 
with a slight transverse polarization at small xp, which turns to longitudinal at 
Xf > 0.6. No Pt dependence of A is observed. Theoretical calculations are needed 
in order to interpret the measurements. 
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CHAPTER 1 

INTRODUCTION 

Since the 1970s the Standard Model has provided a satisfactory description of 
the interactions of all known elementary particles. The underlying theory for de- 
scribing the electromagnetic force in the sub-atomic world is known as Quantum 
Electrodynamics (QED), while Quantum Chromo dynamics (QCD) describes the 
strong force. QED has been tested to be a valid theory by its amazing predictions 
of the lepton magnetic moments and atomic energy spectra. But unlike QED, 
even though QCD was developed following the same fundamental idea of gauge 
invariance and seems to be a straightforward extension of QED, QCD is facing the 
difficulties associated with non-perturbative calculations in the low-energy regime. 
After decades of effort physicists have developed techniques, such as renormaliza- 
tion and resummation, as well as non-perturbative ones, such as effective theories, 
to help solve some of the mathematical difficulties. And nowadays we are able to 
compare many experimental results with QCD predictions, and indeed QCD has 
proved to be the best candidate theory for describing the strong interaction. 

There are many successful examples of QCD. The earhest and most pro- 
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found one is the prediction of the evolution of the structure functions in Deep 
Inelastic Scattering (DIS). Given the quark and gluon distributions at a fixed 
energy-momentum transfer Qq, QCD can actually predict the nucleon struc- 
ture functions at arbitrary using the Evolution Equations fPic 95||. Later in 



collider experiments, Next-to-Leading Order (NLO) QCD calculations predicted 
the inclusive jet-production cross section over several orders of magnitude and 
over a wide range of center-of-mass (COM) energy and jet transverse momen- 



Arn 86 




Ah 91, 


Abe 93 



extracted from the DIS data. Another example is the Drell-Yan process |Pre 70 
The Drell-Yan process, massive lepton-pair production via electroweak quark- 
antiquark annihilation into vector bosons (photon, W^, or Z) and then decay, is 
one of the few processes which have been calculated up to next-to-next-to-leading 
order (NNLO) in perturbative QCD. A recent calculation |[Rij 95|| including some 



of the NNLO terms has shown a good agreement with the data. Another exam- 
ple is the inclusive heavy quark production. Fixed-target studies of heavy-flavor 
production have provided a wealth of data. Total cross sections, single-inclusive 
distributions, correlations between the quark and the antiquark have been mea- 
sured in both hadro- and photoproduction. All experimental results are in quali- 
tative agreement with perturbative QCD calculations. A detailed comparison of 
the fixed-target data and NLO QCD predictions can be found in |p^ri 97| . 



Throughout the entire thesis the following symbols are used to describe the 
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kinematic variables: 

S: center-of-mass energy of the beam-target system. 

m: rest mass of the dimuon pair. 

Pt'- transverse momentum of the dimuon pair. 

xp'- dimensionless longitudinal momentum of the dimuon pair. It is defined 
as the pair longitudinal momentum divided by its maximum kinematically 
allowed value PL,max in the beam-target COM frame. It relates to the Bjorken x, 
the fraction of the hadron momentum carried by the parton in the hadron boosted 
to the infinite-momentum frame, of the beam parton xi and of the target parton 
X2 by Xpi^ — m'^/S) = Xi — X2. 

9 and 0: polar and azimuthal angles of the dimuon pair; described in section 

1.2.3. 

1.1 Failure of Perturbative QCD in Charmonium Produc- 
tion 

While QCD has provided successful descriptions of many aspects of the experi- 
mental data, there are still some phenomena which could not be described. The 
production of charmonium at large transverse momentum is one of such processes: 
the observed production cross sections have differed from QCD predictions by 
more than an order of magnitude, even though in | Fri 97 ] it has been shown that 



the charm production total cross section can be calculated from QCD. This dis- 
crepancy between experiment and theory has revealed a more complicated picture 



for heavy quarkonium production. 

Quarkonium production was conventionally calculated based on the color- 
singlet model (CSM) before 1993. However this model has failed to describe 
charmonium-hadroproduction data |[Sch 94||. In hadropro duct ion of charmonium 



at fixed target energy, \^ < 50 GeV, the ratio of the number of J/\l''s pro- 
duced directly to those arising from decays of higher charmonium states is under- 



predicted by at least a factor of five [Van 95|. At Tevatron collider energies, the 



excess of direct ip' production compared to the CSM prediction is a factor of 30 



Bra 94|, [Roy 94|| . This excess has been referred to as the ifj' -anomaly. Figure ITTT 



compares the CDF ^jj' data and some theoretical predictions. Figure [L^ shows 
the comparison with the fixed-target data. 

1.2 Developments of Theoretical Models 

Data from the Tevatron have revealed that the production rate of tp' at large 
transverse momentum is more than an order of magnitude larger than the early 
theoretical predictions. These results can be understood by taking into account 
two more mechanisms. The first is the realization that fragmentation must domi- 
nate at large transverse momentum, which implies that most charmonium in the 
large region is produced by the hadronization of individual high-p^ partons. 
The second is the development of a factorization formalism for quarkonium pro- 
duction based on non-relativistic QCD (NRQCD) that allows the formation of 
charmonium from color-octet cc pairs to be treated systematically. In this section 




5 10 15 20 

Pt (GeV) 

Figure 1.1: CDF data on the differential cross section for prompt i/j's that do 
not come from Xc decays as a function of pt- The curves are the LO predictions 
of the color-singlet model (dashed curve), predictions including fragmentation in 
the color-singlet model (dotted curve), and including contributions from gluon 
fragmentation via the color-octet mechanism (solid curve) with the normalization 
adjusted to fit the CDF data. (Taken from Pra 96|| ) 
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Figure 1.2: J/'i/' production cross section {xp > 0) in proton-nucleon collisions. 
Solid line: Fit to the total cross section including radiative feed-down from the 
Xc states and ijj'. Dashed line: Direct J/ip cross section. Dotted line: Direct J/ip 
production, color-singlet only. (Taken from Pnk 97|| ) 
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we will summarize these theoretical developments. A more complete review of the 
development was given by Braaten et al. ||Bra 96 1. 

1.2.1 Color Singlet Model 

A thorough review of the applications of the color-singlet model to heavy- quarkonium 
production was given by Schuler [^ch 94|| . To describe the color-singlet model, we 



can think of the production of charmonium as proceeding in two steps. The first 
step is the production of a cc pair, and the second step is the binding of the cc 
pair into a charmonium state. 

We first consider the production of the cc pair. The cc pair must be produced 
with relative momentum that is small compared to the mass of the charm quark 
in order to have a significant probability to be bound together. Assuming that the 
c and c do not exist in the initial state, the production of a cc pair must involve 
virtual particles which are off-shell by amounts of order rric or larger. This part 
of the amplitude is called the short-distance part, because the spatial separation 
of the c and c is of order l/rric or smaller. On the other side, the formation of 
the bound state is considered to be the long-distance part of the amplitude. The 
total amplitude of charmonium production is expected to be dependent on the 
charmonium state H and on the quantum numbers of the cc bound pair. 

For any charmonium state, the dominant Fock state is a color-singlet cc pair 
in a definite angular-momentum state. We introduce the following notation, for 



example, the dominant Fock state for the J/^p is \cc{l, Si)), while for the 



cj 



it is |cc(l,^Pj)). The color states are denoted by 1 for color-singlet and 8 for 
color-octet, and the angular momentum states are denoted using the standard 
spectroscopic notation ^'^+^Lj. The color-singlet model requires that, only the 
cc pair in a color-singlet ^^~^^Lj state can bind to form the charmonium with 
|cc(l^'^+^Lj)) as the final Fock state. 

The color-singlet model has enormous predictive power. The cross section 
for producing a quarkonium state is predicted in terms of a single nonperturba- 
tive parameter for each orbital angular momentum multiplet. The amplitude for 
producing a color-singlet cc pair with small relative momentum (the short dis- 
tance part) can be calculated using perturbative QCD, while the long-distance 
part parameters can be determined from experiments and are expected to be 
process- independent. Thus the long-distance parameters determined from decays 
of the charmonium states can be used to predict the normalized production rate 
of charmonium states. 

We should keep in mind that the color-singlet model is only a model. The 
most basic assumption, the factorization picture, has never been proven to be 
correct, and the relativistic corrections which account for the relative velocity of 
the quark and antiquark are neglected. The color-singlet model also assumes that 
a cc pair produced in a color-octet state will never form the final charmonium. 
However it might be possible that a color-octet cc pair can transit to a color-singlet 
state by radiating soft gluons. We will include the color-octet mechanism in the 
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coming section. 

1.2.2 Gluon Fragmentation 

The first major conceptual advance in recent theoretical developments of quarko- 
nium production was the idea of "fragmentation." Fragmentation is the formation 
of a hadron within a jet produced by a parton with large transverse momentum. 
But here this term is used to include general hadronization processes. 

The real revolution about the fragmentation mechanism is the realization 
that a colored parton, generally a gluon, can result in a color-singlet final state via 
soft-gluon emissions. This possibility was not considered in the conventional wis- 
dom. Once it was accepted, the color-octet cc state could also result in the color- 
singlet final quarkonium by the same argument, and thus the contributions from 
color-octet components become possible, as opposed to the color-singlet model. 

When the GSM includes the contributions from gluon-fragmentation, its pre- 
diction qualitatively agrees with the shape of the CDF ip' data, but is still off in 
normalization by an order of magnitude, in the ip' px differential cross section. 
The prediction from the CSM failed completely in the high region without 
including the gluon-fragmentation mechanism [ |Bra 95[ |. Figure |0| shows the pre- 
dictions and data. 
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1.2.3 Color-Octet Mechanism 

The second major conceptual advance is to realize that the color-octet mechanism 
can be important. Contrary to the basic assumption of the color-singlet model, a 
cc pair that is produced in a color-octet state can bind to form the charmonium 
final state. 

By including the contribution from the color-octet object in the matrix el- 
ement, one can make the prediction agree well with the experimental data by 
leaving the relative size of the color-octet contribution as an adjustable param- 
eter. Its verification now requires considering quarkonium production in other 
processes in order to demonstrate process-independence of the long-distance part 
of the color-octet matrix element. Now the data available from different processes 
are CDF data, fixed-target data, and photo-production data. The size of color- 
octet contributions from these data are not obviously in agreement with each other 
and more sophisticated explanations are needed. 

Other problems associated with the color-octet mechanism are the discrepan- 
cies with the Xci/Xc2 production ratio and the J/ip{ip') polarization: the Xci/Xc2 
production ratio remains almost an order of magnitude too low, and the predicted 
transverse polarization of the J/ ip and ip' is too large compared to the existing pion 
data in fixed-target experiments ||Ben 96||. All of these suggest that higher-twist 



effects may be substantial even after including the octet mechanism. 

A polarization measurement is a crucial test for the color-octet mechanism. 
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(a) Color-Singlet Feynman Diagrams 




(c) Color-Octet Feynman Diagrams 



Figure 1.3: Feynman diagrams of (a) Color-Singlet mechanism, (b) Gluon frag- 
mentation, (c) Color-Octet mechanism. 
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Since the octet production matrix elements of NRQCD lead to a polarization 
pattern different from the CSM, a polarization measurement can provide us with 
significant information on quarkonium production. For example, J/ip's produced 
at the Tevatron at large pt are predicted to be almost fully transversely polarized, 
i.e. X{J/tp) ~ 1 ||Cho 95||, as a result of production via gluon fragmentation. At 



smaller pt, the J/ip's are predicted to be produced essentially unpolarized around 



Pt ~ 5 Gev |Pen 97|| . The observation of this polarization pattern would test 
the underlying theory (the Factorization Approach). To limit the introduction, 
we will concentrate our attention on fixed-target experiments from now on since 
FNAL E866 is a fixed-target experiment. 

The polarization of the quarkonium, measured by analyzing the angular dis- 
tribution of the quarkonium decay products in its rest frame, is of the form 

dcr/rfcos6' ~ 1 + Acos^6l (1.1) 

where 6 is the polar angle measured in the rest frame of the quarkonium. The 
quarkonium rest frame is well specified except for arbitrary three-dimensional 
rotations. The CoUins-Soper frame \ \Joi 77|| , in which the Z-axis is defined to 



be parallel to the bisector of the angle between the directions of the interacting 
hadrons in the quarkonium rest frame, is used in this analysis. In all other ear- 
lier fixed-target experiments the Gottfried- Jackson frame, in which the Z-axis is 
defined to be parallel to the incoming beam axis in the quarkonium rest frame, 
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was used. These two frames are equivalent if the quarkonium has zero px [ Kal 86 
For the pt range of fixed-target experiments, of the order of 1 GeV, compared to 
hundreds of GeV of longitudinal momentum, the two frames are approximately 
the same IIGee 98 . 



1.3 Fixed- Target Polarization Experiments and Predictions 

Polarization measurements have been performed for J/ip and ip' production in 
pion and proton scattering fixed-target experiments. From a theoretical point 
of view, the ip' decay has been more extensively studied because all the ip' data 
samples are direct ip's. The observed value of A for ip' is 0.02 ± 0.14, measured at 
^/s = 21.8 GeV in the region xp > 0.25 by Heinrich et al. |Pei 9 If . When studying 
the polarization of the J/ ip decay one has to take the polarization inherited from 
decays of the higher charmonium states Xcj and ip' into account and this leaves 
some ambiguity in the interpretation of the results. In the following sections we 
will only review J/ tp polarization experiments to compare with the E866 results. 

1.3.1 Model Predictions of Polarization at Fixed- Target 
Energies 

1.3.1.1 Color-Singlet Model 

The polarization of J/ip has been calculated from perturbative QCD by Vanttinen 
et al. IJVan 95||. The parameter A in Equation 1.1 was calculated from the cc 



production amplitude and the electric dipole approximation of radiative x decays. 



Figure |1.4| a shows the predicted values of the parameter A in Equation 1.1 
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Figure 1.4: Leading-twist predictions of the parameter A in the decay angular 
distribution of J/ip's produced in pion-nucleon colhsions at beam energy Eiab = 
300 GeV, plotted as a function of xp. (a) The three solid curves show the decay 
distributions of J ftp's produced via radiative decays of the X2 and Xi states and 
"directly" in gluon fusion. The dashed curves show the effect of smearing the 
transverse momentum distribution of the beam parton by a Gaussian function 
exp[-(A;j_/500 MeV)^]. (b)The combined decay distribution of all J/V^'s, including 
contributions from xi,2 decays and direct production, is shown here. The lower 
curve shows the effect of adjusting the relative normalization of the different con- 
tributions to their measured values by appropriate K-factors. The dashed curve 
shows the effect of transverse-momentum smearing and K-factors adjustments. 
The data are from the Chicago-Iowa-Princeton (full circles) and E537 (open cir- 



cles) experiments. (Taken from [[Van 95 
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in the Gottfried- Jackson frame as a function of xp, for the direct J/ip and the 
Xi,2 — J/ip + 'J processes separately. The dashed hnes indicate the effect of a 
Gaussian smearing in the transverse momentum of the beam partons. The over- 
all \{xf) including direct and indirect J/ip processes is shown in Figure p,.4|b 



and compared with the Chicago-Iowa-Princeton ||Bii 87|| and E537 data ||Ake 93 



The QCD calculation gives A ~ 0.5 for xp < 0.6, significantly larger than the 



measured value. The lower curve in Figure |1.4| b shows the effect of multiplying 
the partial J/ ip cross section with the K-factors obtained from experiments. The 
discrepancies between the calculated and measured values of A suggest that the 
leading-twist processes considered in the calculation are not adequate for explain- 
ing charmonium production. 

1.3.1.2 Color-Evaporation Model 



The color-evaporation model [[bri 77| , [Hal 77| | assumes that the cc pair in ^5*1 state 



can transit to ^Sq state via soft gluon emission, so J/ip is always produced un- 
polarized. In this model the color and spin quantum numbers of the cc pair are 
irrelevant. The fraction of the cc pairs bound into J/ip is described by a phe- 
nomenological parameter fj/^. 

The color-evaporation model is considered to be an over-simplified model, 
because it is not concerned with the details of the particles which initiate the 
reaction. The evident failure is the prediction of the fraction of J / ip coming from 
Xc decays. According to the color-evaporation model, the fraction of J/ip coming 
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from Xc decays should be process-independent. But the experimental data both 
in fixed-target experiments in pN and vrN collisions and also in pp collisions at the 
Tevatron gather around a central value of 0.3-0.4, while in 7-p collisions an upper 
limit of 0.08 was obtained IIBar 871 . 



Since this model gives trivial prediction on J/ if) polarization and fails in 
predicting ratios of quarkonium production, we will not discuss this model in 
later discussions. 

1.3.1.3 Non-Relativistic QCD 

The polarization of J/ip has been calculated in non-relativistic QCD by Beneke 
and Rothstein ||Ben 96 1. The production cross section for a quarkonium state H 



in the process 

A + B — >H + X (1.2) 

can be written as 



(Th = Y1 j dxidx2fi/A{xi)fj/B{x2)a{ij H) (1.3) 
a(.j-i7) = EC^|Q[„](0f) (1.4) 

n 

In Equation 1.3 the summation sums up the contributions by all partons 
in the colliding hadrons, and the fi/A and fjjB are the corresponding parton 
distribution functions (PDF). The coefficients C'^gj^] in Equation 1.4 describe the 
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production of a quark-antiquark pair in a state n and have expansions in asi^mq). 
The parameters (O^) describe the subsequent hadronization of the QQ pair into 
the quarkonium state H. It is important to test the universahty of the production 
matrix elements (O^) because this is an essential prediction of the factorization 
formula (1.4). 

In the calculation of Beneke and Rothstein, the following intermediate cc 
states are considered: (L^S*!), (8/ 6*0), (8,^Pj), and (8,"^ 5*1). For each intermediate 
state the ratios of longitudinal to transverse polarized quarkonia were computed. 
To obtain the total polarization, the various subprocesses have to be weighted 
by their partial cross sections. Weighting all subprocesses by their partial cross 
sections and neglecting the small ip' feed-down, a sizable polarization is obtained: 

0.31 < A < 0.63 

However the existing data show no sign of polarization. Thus NRQCD in- 
cluding the color-octet contributions also gives a wrong prediction on the J/ip 
polarization problem, and one has to seek for explanations from higher-twist pro- 
cesses. 

1.3.2 Fixed-Target J/i/j Polarization Experiments 
1.3.2.1 E537 

Fermilab experiment E537 has measured the differential cross section da/dcos6 

for J/ijj production in 7r~N interactions and in pN interactions at ^/S = 15.3 

GeV in the region xp > [ Ake 93 |. Fitting the angular distribution to the form 
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Figure 1.5: cos6' distribution for E537 pW data. (Taken from ||Ake 93|| ) 



of Equation 1.1, A = -0.115 ± 0.061 for p and A = 0.028 ± 0.004 for tt" were 
obtained. The data sample used to obtain this resuh contained 12530 J/ip events 
produced by the p beam and 33820 J/ip events by the 7t~ beam. 
1.3.2.2 E672/E706 

Fermilab experiments E672/E706 have measured the differential cross section 
da/dcos6 for J/^p production in 7r~Be collisions at = 31.5 GeV in the re- 
gion 0.1 < xp < 0.8 ||Gri 96||. Fitting the angular distribution to the form of 



Equation 1.1, A = —0.01 ± 0.08 was obtained. The data sample used to obtain 
this result contained 9600 J/ip^s. 
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Figure 1.6: cos6' distribution for E537 vr W data. (Taken from |[Ake 93|| ) 



1.3.2.3 E771 



Fermilab experiment E771 has measured the differential cross section da/dcos6 
for J/ijj production in pSi coUisions at \/S = 38.8 GeV in the region —0.05 < 



xp < 0.25 ||Ale 9711 . This is the only published polarization measurement for J/ijj 
produced with a proton beam. Fitting the angular distribution to the form of 
Equation 1.1, A = —0.09 ± 0.12 was obtained. The data sample used to obtain 
this result contained 11660 J/ip^s. 
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Figure 1.7: xp, Pt, and /cos6 distribution for the E672/706 vr Be data. (Taken 
from ||Gri 96|| ) 
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Figure 1.8: cos 6' distribution for E771 pSi data. (Taken from [|Ale 97|| 



1.3.2.4 Chicago-Iowa-Princeton 



A dedicated J/ip decay angular distribution measurement was performed at Fer- 



milab using a 252 GeV pion beam incident on a tungsten target ||Bii 87|| . The 
data sample contains 1600000 J/ip events from a vr" beam and 600000 J /ip events 
from a vr^ beam. The data are in the kinematic range xp > 0.25 and Pt < 5.0 
GeV. To determine the J/ ip decay angular distribution, the data were divided into 
fifteen regions of xp-, five regions of cos^, and five regions of in the kinematic 
range xp > 0.25, —1 < cos6' < 1, and — vr < (p < n. For each bin of xp, cos 9, 
and (f) the raw fj^^fi~ mass distribution was fitted by a seven-parameter function 
involving a Gaussian distribution for the J/ ip and tp' and a quadratic polynomial 
plus an exponential of a first-order polynomial for the continuum background. 
The number of J/ip^s in the 375 bins of xp, cos6, and (p were then corrected for 
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Figure 1.9: The xp dependence of the parameters fit to the J/tp decay angular 
distribution in Equation 1.5. (Taken from |[Bii 87|| ) 



acceptance and for each of the fifteen regions of xp the J/ip angular distribution 
was fitted by the general form | Lam 78 1 



d'^a/d cos 9d(j) ~ 1 + A cos^ 6 + fi sin 26' cos(f) + -u sin^ 6 cos '. 



:i.5) 



The A, fi, and u are consistent with zero over a wide range of xp- Note that A 
approaches —1 at high xp- This behavior was also observed in a Drell-Yan con- 



tinuum production experiment [|Ale 86|| . Table [L^ summarizes the experimental 
and theoretical results of J/ip polarization. 
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Table 1.1: Summary of the experimental and theoretical results. 



Experiment reaction xp range A 

E537 p + W 15.3 GeV xp > -0.115 ±0.061 

E537 TT- + W 15.3 GeV xp > 0.028 ± 0.004 

E672/706 TT- + Be 31.5 GeV 0.1 < xp < 0.8 -0.01 ± 0.08 

E771 p + Si 38.8 GeV -0.05 < xp < 0.25 -0.09 ± 0.12 

CIP TT + W 21.7 GeV 0.25 <xp<l.Q ^0, ^ -1 at large xp 



Theory xp range A 



GSM 


Xp 


> 


~ 0.25 


GEM 


Xp 


> 





NRQGD 


Xp 


> 


0.31 < A < 0.63 



1.4 Fermilab E866 Measurement 

E866 at Fermilab was designed to measure the u/ d asymmetry in the nucleon sea. 
After the run ended in March, 1997, additional measurements were performed in 
the run extension period. Two major topics in the run extension were angular 
distribution of the J/ip decay and nuclear dependence of J/ip production. The 
work presented here is based on the data sample collected during a four week 
dedicated beam-dump run, from which the angular distribution of the J/ip decay 
in the dilepton channel was studied. This angular-distribution measurement is 
unique since no high-statistics proton- induced data exists. Also, the J/ip produc- 
tion diagrams are different for pN and ttN interactions. A total of 10 million J/tp^s 
(with ~ 1% of unseparated ■0"s) were collected, and the kinematic coverage of the 
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data extends over xp > 0.2, px < 5GeV, and —0.95 < cos^ < 0.95. The quantity 
of the data sample has allowed us to present the A parameter in Equation 1.1 in 
seven regions oi xp and four regions of pr- The results could provide a test of 
the color-octet mechanism, and hopefully will improve our understanding of the 
higher-twist effects. 
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CHAPTER 2 



EXPERIMENTAL APPARATUS 



The experiment E866 was performed at the Meson-East experimental area of 



Fermi National Accelerator Laboratory. The spectrometer, shown in Figure ^^1 
was a modified version of the E605/E772/E789 spectrometer. This spectrometer 
was designed to detect dimuon events with forward xp, though certain combi- 
nation of target position and analyzing magnet settings allows finite negative- 
xp acceptance. The spectrometer primarily consisted of three dipole magnets, 
seven hodoscope planes, eighteen drift-chamber planes, and three proportional- 
tube planes. The hodoscope planes were used to provide the trigger information, 
the drift-chamber planes were used to find the trajectories, and the proportional- 
tube planes, which were also part of the trigger system, were used to identify 
muons. The SMS magnet measured the momentum of the muon pairs while the 
SMO and SM12 magnets allowed us to select the desired mass range. The charged 
particles produced in the target were split according to the sign of their charges 
while going through the set of three magnets. 

A thick hadron absorber wall was installed in SM12 for the experiment so 
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Figure 2.1: The E866 spectrometer. 

that the long-hved hadrons (mainly pions) can be stopped before hitting station 
1. The absorber wall consisted of Cu, C, and (CH2)n blocks, and gave a hadron 
attenuation factor of e~^°. Not shown in Figure 2.1 is the copper beam dump 
sitting in front of the hadron absorber wall. The beam dump was used to stop 
the 800 GeV proton beam. For the angular distribution measurement we used 
the beam dump as the target. There were also an electromagnetic and a hadronic 
calorimeter and a ring-imaging Cherenkov counter (RICH), which, however, were 
not operating during this experiment. The RICH counter was filled with helium 
bags to reduce multiple scattering. This chapter will discuss only these compo- 
nents that were needed for this study. 

Throughout, we will make reference to the spectrometer-fixed coordinate 
system. The E866 coordinate system aligns the Z-axis horizontally with the ac- 
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celerator proton beamline and the Y-axis with the vertical direction. The X-axis 
is then chosen to form a right-handed Cartesian system. The positive Z direction 
is chosen to be same as the beamhne direction, which is also referred to as "down- 
stream," and the positive Y direction is chosen to be up. The origin is located at 
the center of the upstream face of the SM12 spectrometer magnet. 

2.1 Accelerator and Beam 

The high-energy proton beam was produced in the Tevatron, which is a supercon- 
ducting proton synchrotron. Protons were first accelerated by a pre-accelerator 
up to about 700 KeV. These protons were then accelerated in a linear accelerator 
to about 400 MeV. Subsequently, a Booster Ring boosted the proton energy to 8 
GeV. Protons were then injected into the main ring, located in the same tunnel as 
the Tevatron but constructed from conventional magnets, in which protons could 
reach an energy of 400 GeV before being transferred into the superconducting 
ring, where the protons were accelerated to 800 GeV. After the proton beam was 
accelerated to 800 GeV, it was extracted and split by the switchyard for sending 
three streams of proton beams to the Meson, Neutrino, and Proton beam lines 
for the fixed-target experiments. 

Protons in a spill were bunched into RF buckets separated from one another 
by 18.9 ns, with bucket length ~ 1 ns. Each spill contained about 10^ buckets. 
This small scale beam structure was due to the Tevatron accelerating radio fre- 
quency of 53 MHz. A square wave signal at this frequency, called the RF clock, 
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was used to synchronize the E866 electronics with the Tevatron beam structure. 
It took ~ 25 seconds to accelerate a fill of protons up to 800 GeV. These protons 
were then extracted from the accelerator for 20 seconds. After that, the super- 
conducting magnets ramped down for 15 seconds. The entire cycle time was 
approximately one spill per minute. Typical proton intensities in the Tevatron 
were 1 — 2 * 10^^ protons per spill. 

Within the Meson three-way split divided the proton beam between 

the Meson- East line and the rest of the Meson lines. To monitor the beam in- 
tensity, luminosity, position, and beam-spot size, several beamline detectors and 
monitors were used. During the beam-dump running mode the typical beam in- 
tensity was 6* 10^° protons per spill. The beam intensity was monitored by an ion 
chamber located in the ME3 sector (IC3), a secondary-emission monitor located 
in the MEG sector (SEM6), and a beam Cherenkov monitor. Both the size and the 
position of the beam were monitored by segmented wire ion chambers (SWICs) 
and the Beam- Position Monitor (BPM). The beam luminosity was monitored by 
the AMON and WMON scintillation counters, which were installed at about 85 
degrees from the target position. 

2.2 Beam-Dump Target 

For this study the beam dump itself was the target. The dump was suspended 
from two of the central magnet inserts inside the SM12 magnet, beginning at Z 
= 68 inches, extending 168 inches downstream, and ending at Z — 236 inches. In 
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Figure |2.2| a picture of the beam dump and absorber wall is shown. The beam 
dump was made of pure copper with cooling water tubes running through the 
sides. There was a 12-inch deep rectangular-shaped hole in the center of the 
upstream face of the dump to help contain backscattered particles, so actually 
the beam protons did not hit the dump until Z = 80 inches. This still left 156 
inches thickness of copper, which is equivalent to 26.5 interaction lengths for p + 
Cu collisions, to stop the protons and secondary particles. The probability for a 
primary proton to punch through the entire dump was less than 4 * 10^^^. Most 
primary protons would interact within the first few interaction lengths. Secondary 
particles from the primary interaction would further interact to form showers and 
eventually be stopped in the dump, but the high-energy muons produced would 
penetrate the entire dump with little interaction since muons are not strongly- 
interacting particles. However, these muons would still lose energy and suffer 
multiple scattering on their way through the dump, and thus added uncertainties 
to the reconstructed xp and px- 

2.3 Spectrometer Magnets and Absorber Wall 

Two dipole magnets, SM12 and SMS were used in this study (the current of SMO 
was set to zero during the beam-dump data taking). The magnetic fields of these 
magnets were oriented horizontally. The field strengths of the dipole magnets 
could be configured to optimize the mass acceptance for J/ip^s. The bending 
magnet used for this study was SM12. The length of the SM12 magnet, made of 
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Elevation View (Y) 

Figure 2.2: The SM12 acceptance defining magnet with the absorber wall is shown. 
The magnet coils and iron return yoke are only partially shown. The beam dump 
is also not shown in the plan view. 

iron, was 567 inches. The magnet produced an average horizontal field of up to 
1.3 Tesla at a maximum current of 4000 Amperes. This corresponds to a 7-GeV 
transverse momentum kick to the charged particles which traveled through its 
entire length. In this study, SM12 was set to 2800 Amperes and 2040 Amperes 
during two separate data-taking periods, delivering a transverse-momentum kick 
of 4.2 GeV and 3.1 GeV, respectively. 

The momenta of the muon tracks were measured by the analyzing magnet 
SM3. The location of the SM3 magnet was between Station 1 and Station 2, as 
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shown in Figure SMS delivered a transverse momentum kick of 0.914 GeV 
to the charged particles traveling through when operated at its maximum current 
of 4260 Amperes. The field was uniform enough so that the reconstruction of a 
particle trajectory through the field volume can be described by a single bend- 
plane approximation. 

The absorber wall was located inside the SM12 magnet directly behind the 
beam dump. It filled the SM12 magnet completely in the x and y direction. 
The absorber wall was constructed of one 24-inch section of copper, three 27-inch 
sections of carbon graphite, one 27-inch section of carbon-polyethelene compound, 
and two 36-inch sections of polyethelene. 

Both magnets were filled with helium bags to minimize the multiple scatter- 
ings of the muons. 

2.4 Detector Stations 

There were four detector stations in the E866 spectrometer, denoted as Station 1 
to 4. Station 1-3 each consisted of hodoscopes and drift chambers, while Station 4 
consisted of hodoscopes and proportional tubes. Those stations record the passage 
of charged particles in space and time across their active area. Together with the 
information provided by the magnet field maps, this allowed the 4-momentum of 
the individual tracks to be reconstructed. Stations 1-3 were used for triggering 
and tracking, while Station 4 was used for muon identification and triggering. 
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2.4.1 Drift Chambers 

Each one of Stations 1-3 consisted of 6 planes of drift chambers. The 6 planes 
were arranged in pairs with parallel wire orientation. The second plane of a pair 
had its wires offset by half the cell size of the drift chamber. The upstream 
plane of each pair was denoted as the "unprimed" plane, while the downstream 
plane was denoted as the "primed" plane. The Y-Y' pair of each station held the 
wires horizontally to measure the Y-intercept of the tracks, while the V-V and 
U-U' chambers had their wires tilted at —14 (a slope of —0.25) degrees and +14 
(a slope of 0.25) degrees from the X-axis respectively. These planes determined 
the X-intercept of the track and also provided a check on the Y-intercept. The 
configuration of the drift chambers are given in Table [2?!. 



The drift chambers were all operated with a gas mixture of 49.7% argon, 
49.6% ethane, and 0.7% ethanol by volume, which was mixed at a constant 25 °F. 
The Station-1 anode wires were made of gold-plated tungsten wire, while Stations 
2 and 3 used silver-coated beryllium-copper wires. All the anode wires were 25 fan 
in diameter. The cathode wires for all three stations were silver-coated beryllium- 
copper wire with a diameter of 62.5 /im. The absolute operating voltages were 
between 1700 and 2200 volts, which gave a typical drift velocity about 50 /im/ns. 

The signals of these chambers were read out by a fast amplifier and discrim- 
inator system. Single-hit TDCs (Time-to-Digital Converters), which only record 
the first hit on the wire during an event, were used to measure the drift time. The 
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Table 2.1: Drift chamber parameters. The unit length is one inch. 



detector 


Z-position 


No. of wires 


cell size 


aperture (XxY) 


operating voltage 


VI 


724.69 


200 


0.25 


48x40 


+1700 


VI' 


724.94 


200 


0.25 


48x40 


+1700 


Yl 


740.81 


160 


0.25 


48x40 


+ 1700 


Yl' 


741.06 


160 


0.25 


48x40 


+ 1700 


Ul 


755.48 


200 


0.25 


48x40 


+1700 


Ul' 


755.73 


200 


0.25 


48x40 


+1700 


V2 


1083.40 


160 


0.388 


66x51.2 


-2000 


V2' 


1085.52 


160 


0.388 


66x51.2 


-2000 


Y2 


1093.21 


128 


0.40 


66x51.2 


-2000 


Y2' 


1095.33 


128 


0.40 


66x51.2 


-2000 


U2 


1103.25 


160 


0.388 


66x51.2 


-1950 


U2' 


1105.37 


160 


0.388 


66x51.2 


-1975 


V3 


1790.09 


144 


0.796 


106x95.5 


-2200 


V3' 


1792.84 


144 


0.796 


106x95.5 


-2150 


Y3 


1800.20 


112 


0.82 


106x91.8 


-2200 


Y3' 


1802.95 


112 


0.82 


106x91.8 


-2200 


U3 


1810.24 


144 


0.796 


106x95.5 


-2200 


U3' 


1812.99 


144 


0.796 


106x95.5 


-2200 
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combination of good hits together with their associated drift times in all three 
views gave a "triplet" hit for a station. The bank of the triplets was saved to 
provide information for the track reconstruction. 

2.4.2 Hodoscopes 

Associated with the drift-chamber planes, there were also hodoscope planes in 
each tracking station. These hodoscopes provided fast tracking signals for use 
in triggering. In Stations 1, 3, and 4 there were two hodoscopes planes which 
measured the X and Y intercepts of the tracks, while in Station 2 there was only 
one hodoscope plane. Each hodoscope plane was arranged into two half-planes 
of parallel scintillator paddles, which were attached to photomultiplier tubes via 
plexiglass light guides. During operation, each paddle only gave a single bit of 
signal (one or zero). 

The hodoscope planes were named according to the tracking station they 
belonged to, preceeded by X or Y depending on the orientation of the paddles. 
For example, 'Y3 hodoscope" referred to the Station-3 hodoscope plane in which 
2x13 scintillator detectors was positioned horizontally and separated into left and 
right side. The parameters of the seven hodoscope planes are given in Table P^ . 

2.4.3 Proportional Tubes 

Station 4 was also called the muon station. It was located downstream of the 
calorimeters and consisted of two hodoscope planes (Y4,X4) and three propor- 



34 



Table 2.2: Hodoscope plane layout. Dimensions are in inches. 



detector Z-position No. of counters cell width aperture XxY 



Yl 


769.78 


2x16 


2.5 


47.50x40.75 


XI 


770.72 


12x2 


4.0 


47.53x40.78 


Y2 


1114.94 


2x16 


3.0 


64.625x48.625 


X3 


1822.00 


12x2 


8.68 


105.18x92.00 


Y3 


1832.00 


2x13 


7.5 


104.00x92.00 


Y4 


2035.50 


2x14 


8.0 


116.00x100.00 


X4 


2131.12 


16x2 


7.125 


126.00x114.00 



tional tube planes(PTYl, PTX, PTY2). Each of the three proportional tube 
planes had two layers of 1 x 1-inch cells. These two layers were offset by a half-cell 
spacing to cover the dead region between the adjacent cells. The proportional 
tubes used the same argon/ethane/ethanol gas mixture as the drift chambers. 
To minimize the probability of hadron punch-through, an absorber wall (3 feet of 
zinc and 4 inches of lead) was placed between the calorimeter and the muon detec- 
tor. Furthermore, 3-foot thick concrete walls were placed between PTYl and X4, 
and between PTX and PTY2. This provided a total of 16.6 interaction lengths 
upstream of Y4. Thus the only charged particles which could reach Station 4 
detectors were the muons. Signals from the cell of the proportional tubes were 
amplified and shaped by the attached pre- amplifier/discriminator cards. Signals 
exceeding the threshold voltage were sent to the Coincidence Registers (CRs) to 
indicate the arrival of muons. The parameters of all proportional tube planes are 
given in Table 1^ . 



35 



Table 2.3: Proportional tube parameters. All dimensions are in inches. 



detector 


Z-position 


No. of wires 


cell size 


aperture XxY 


PTYl 


2041.75 


120 


1.0 


117x120 


PTX 


2135.875 


135 


1.0 


135.4x121.5 


PTY2 


2200.75 


143 


1.0 


141.5x143 



2.5 Trigger System 

2.5.1 Trigger-System Hardware 

A new trigger system was implemented for E866 data taking [ [Haw 98|| . A block 
diagram of most of the trigger system is shown in Fig. 2.3. 
2.5.1.1 Trigger Matrix Module 

Scintillator counters were used to provide input signals for the trigger system. 
Signals from the photomultiplier tubes attached to the scintillator counters were 
brought to LeCroy 4416 16-channel discriminators via coaxial cables. Each dis- 
criminator output was synchronized to the accelerator RF signal and shaped to 
a fixed width of 15 ns by pulse stretchers, and then fanned out to Coincidence 
Registers (CRs), Terminator/OR modules, and the Trigger Matrix (TM) modules. 

The Trigger Matrix modules were the core of the trigger system. The pulse- 
stretcher outputs were grouped as a half-bank (right or left) of Yl, Y2, and Y4 
hodoscope planes. These groups were sent as inputs to the Trigger Matrix to 
identify tracks of muons originating in the target. Only Y-view hodo-roads were 
used here, because events with different Z positions and momentum had different 
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the spectrometer are shown. 



"roads" in the Y view under the deflection of the bending magnets. This Trigger 
Matrix was conceptually a lookup table loaded to a set of six 256 x 4-bit ECL 
SRAM chips. All tracks of interest defined a set of valid "roads" going through 
the hodoscope planes. These roads were identified using a Monte Carlo simulation 
for muon tracks coming from the target, and thus a "map" of these roads was 
produced. The "map" was then written into a disk file, which was loaded into 
the SRAMs by the Data Acquisition (DAQ) online program during the start-run 
stage. While taking data, the hits on Y2 and Y4 were combined to form an 
"address" whose content were the predictions on Yl, which were compared to the 
actual Yl hit pattern. Any coincidence in the comparison generated the Trigger 
Matrix output. 

There were four sets of Trigger Matrix modules called MUL, MUR, MDL, 
and MDR. They covered different types of valid muon tracks, namely, up-left, up- 
right, down- left, and down- right, respectively. For finding the target muon pairs, 
the coincident combination of an up and a down track was required. The output 
signals of the Trigger Matrix modules were then sent to the Track Correlator for 
further triggering determination. 
2.5.1.2 Track Correlator 

The Track Correlator (TC) modules were designed by Texas A&M University 
Gag 98||. These programmable modules were used to filter specific combinations 



of Trigger Matrix, Terminator /OR, and S4XY [42] outputs to trigger on an event. 
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However, during beam-dump data taking only output signals from the TMs were 
of interest. Four 16-bit patterns, according to the desired trigger conditions, 
were preloaded into a 2^^ x 4 bit SRAM chip inside a TC during the start-run 
stage. Whenever the output combination of the TM modules matched one of 
the preloaded bit patterns, the TC would, prescaled to the desired frequency and 
synchronized with the RF clock, send out a signal to the Master Trigger OR 
module to notify the arrival of an interesting event. Each SRAM chip could be 
programmed with up to four independent trigger conditions. 

There were three main Track Correlators which were able to trigger on an 
event to start the DAQ. The first TC, called Physics TC A, was programmed to 
select two-tracks events, hke-sign or unlike-sign. The second TC, named Physics 
TC B, was designated for left-right efficiency studies, which involved the use of 
X hodo planes and single muon events. The third module, which was used to 
trigger on cosmic rays to diagnose the trigger and DAQ systems, was called the 
Diagnostic TC. It also provided the measurements for scintillator efficiencies. 

If the inputs to the Track Correlator fulfilled the triggering criteria, a Trigger 
Generate Input (TGI) signal would then be sent to the Master Trigger OR by the 
TC. The Master Trigger OR would then synchronize this trigger signal and the 
drop of DAQ System Busy with the RF clock to send out a Trigger Generate 
Output (TGO). Triggers were thus inhibited during event readout; the difference 
of TGI and TGO counts would provide information on readout dead time. 
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Table 2.4: Correspondence of the SRAM chip bit to various input sources of the 
Physics TC A,B module 

"signal Origin TC SRAMbiT 



S4UL1 


bitO 


S4UL2 


bitl 


S4DL1 


bit2 


S4DL2 


bits 


S4UR1 


bit4 


S4UR2 


bits 


S4DR1 


bite 


S4DR2 


bit? 


MUL 


bits 


MDL 


bit9 


MUR 


bitlO 


MDR 


bit 11 


X134L 


bit 12 


X134R 


bit 13 


- unused - 


bit 14 


- unused - 


bit 15 



2.5.2 Trigger Firing Criteria 
2.5.2.1 PhysA Trigger 

Prom Table |2.5| the definitions of PhysAl,2,3,4 triggers are self-explanatory. PhysAl 
trigger required that two tracks went through two diagonally-opposite quarters of 
the spectrometer, while the PhysA3 and PhysA4 required that two tracks went 
through the same side, left or right, of the spectrometer with one track going 
up and the other going down. These tracks were identified as unlike-sign muon 
pairs and were treated as possible candidates of target events. The PhysA2 trig- 
ger required that both tracks went up or down, and thus gave like-sign muon 
pairs. This information was especially important for rate-dependence studies in 
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Table 2.5: Prcscalc factors and trigger descriptions for Physics TC A,B module. 
The represents a logical AND and the "+" represents a logical OR. 



TtI 0'0'PT Tli^TTlP 
J- ± ICjCjCl iicuiic; 


■n-rpcpQlpT* fap+nT* 

LJl. KjijKjOjlXjl. i-OAj uUl 


r\ pcpfi 1"^+ inn 


PhysAl 


1 


(MUL*MDR) + (MUR*MDL) 


PhysA2 


1 


(MUL*MUR) + (MDL*MDR) 


PhysA3 


1 


(MUL*MDL) 


PhysA4 


1 


(MUR*MDR) 


PhysBl 


10 


(X134L*X134R) 


PhysB2 


1000 


MUL + MDL + MUR + MDR 


PhysB3 







PhysB4 








extracting cross sections. 
2.5.2.2 PhysB Trigger 

The Physics TC B was used for recording events for studies. In the PhysBl trigger, 
the symbol "X134L(R)" represented a track that went through the left(right) side 
of XI, X3, and X4 hodoscope planes. The signals fed into the TC B were outputs of 
some Terminator/OR modules, whose outputs represented the logical ORs of the 
signals of the X hodo scintillators. The trigger requirement, X134L*X134R, was 
designed to measure the random muon coincidences. The other trigger PhysB2 
only required a single hit on any of the four quarters to fire. It had a prescale 
factor of 1000 and was used to measure the rate of single muons. 

2.6 Data- Acquisition System 

The Data- Acquisition System could be divided into three parts by functionality: 
event readout, data archiving, and online analysis. The first part was based on 



41 



a Nevis Transport system, the second was a VME-based data-transferring and 
controlling system, and the third was built on the SGI workstations. 

2.6.1 Readout System 



The backbone of the E866 readout was a Nevis transport system [|Kap 82|| . All 
detector subsystems ultimately fed data into the Transport. The subsystems 
included Time-to-Digital Converter (TDC) readouts from drift chambers and Co- 
incidence Registers (CRs) from hodoscopes and muon proportional tubes signals. 
Bus arbitration was maintained by a hard-wired daisy chain, with the bus mas- 
tership determined by the Carry signal. This scheme not only prevented multiple 
subsystems from attempting to place data on the Transport simultaneously, but 
also guaranteed that events appeared on the readout bus in a well-defined order. 
The data bus was 16-bit wide, and the system clock was set to 10 MHz. All the 
data fed into the Transport Bus were then transferred to a VME-based archiving 
system ||Car 91 . 



Upon receipt of the TGO signal, the first module in the Transport Bus Carry 
chain, the Event Generator Source (EGS), would raise the System Busy signal to 
inhibit any further triggers and take control of the Carry signal. The EGS then 
put a special "first-word" into the Transport bus to indicate the beginning of a 
new event in the data stream. After a few more words from the EGS, the Carry 
signal was passed to the first branch of the readout subsystem to begin transferring 
event data into the Transport. 
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Upon receipt of the TGO signal, the EGS module also fanned out "START" 
signals to CRs and TDCs to begin digitizing the pulse signals. For each hit in 
the hodo or prop tube the CR would generate one word in the event output, 
containing the scintillator ID or wire number of the muon proportional-tube hit. 
Each event also contained a record of which trigger caused it to be readout via 
the Trigger Bit Latch (TBL). In addition to the CR's, each TDC would begin 
incrementing a Gray-code counter once every 4 ns upon receiving the START 
signal. The incrementing process would be stopped by the amplified signal from 
the drift chamber. Each hit in the drift chamber would also produce a one-word 
output, containing the wire number of the hit and the Gray-code value of the 
TDC timer. This measured the drift time. All the data were transferred into the 
Transport bus in the Carry chain order. 

2.6.2 Data Archiving System 

Events from the Transport bus streamed into the VME through a pair of "ping- 
ponging" triple-ported VME high-speed memory boards by way of a front-panel 
ECL interface. Interrupt-driven software would initiate DMA transfers of packets 
of events from the high-speed memories across the VME bus into a 128-megabyte 
ring buffer. This buffer was continuously being drained across the VME bus 
into a single-board computer by a concurrent task which performed all the data 
formatting. Prom there, formatted packets of events were queued in a small pool 
for distribution to the taping subsystem, where up to four Exabyte 8mm tape 
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drives would record the data. 

Unlike in E789, the communications to the readout system and the run- 
control capabihties were all built on the VME single-board computer in E866. In 
addition, the scaler data were injected into the data pipeline, from the CAMAC 
system, as regular logical records on a spill-by-spill basis. These scalers included 
target, beam, magnet parameters and counts from varieties of trigger conditions. 
A small fraction of event packets were fanned to the UNIX workstations for online 
data sampling and analysis. This provided the capability of online monitoring. 

For the beam-dump running, the average data-taking rate was about 20000 
events per spill. The average event size during the beam-dump run was 192 16-bit 
words. 

2.6.3 Data Monitoring System 

The E866 online database system was based upon the ADAMO hbrary distributed 
by CERN, with a graphical interface package called "PinKy." The database for 
E866 recorded various data streams, including the beginning-of-run (BOR), end- 
of-spill (EOS) scalers, and beamline data (EPICURE). The online monitoring 
tools included 1) "runstatus," a graphical display of certain critical data (mag- 
net settings, beam intensity, luminosity, live-time, and duty-factor calculations) 
updated at each EOS, 2) "scan," a graphical display of scalers refreshed at each 
EOS, 3) "plot," a plotting tool for monitoring any entity stored in the database, 
4) "review," a tool for fetching data for series of runs for plotting or exporting to 
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the CERN Physics Analysis Workstation software package (PAW) ntuple file, in 
which the interested quantities of an event were stored in an array, and 5) "dd," 
a tool to receive and distribute data to the backend. These advanced monitor- 
ing tools provided the capabilities for us to reconstruct and monitor a fraction 
of events online during the data taking. The shift taker could, for example, see 
the mass spectrum, hits and multiplicities on the drift chambers and hodoscope 
planes, and format errors due to transport readout problems while the data were 
being taken. Thus this capability helped us to diagnose the hardware problems 
and improved the quality of the data. 
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CHAPTER 3 



MUON-TRACK RECONSTRUCTION 



The data recorded in 8mm tapes are logical records and have to be decoded 
before using. The decoding was accomplished by the analysis code developed for 
E866. The most important information stored in the data stream is actually the 
space-time marks of the electronic signals traveling through the detector stations. 
These marks in reality present the trajectories of the charged particles. From the 
trajectories the kinematic quantities of the particles can then be determined if the 
mass of the charged particle and the strength of the magnetic field is known. 

In this chapter we will first summarize the data taken in April of 1997 for 
this study, and then describe in detail the method of track reconstruction in this 
experiment. 

3.1 Data-Set Summary 

During beam-dump data taking the typical beam intensity request was 6E10, 
6 * 10^° protons per spill. The average triggering rate was about 20000 triggers 
per spill and the average event size was 192 16-bit words, as mentioned in the last 
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Table 3.1: Magnet currents and trigger matrix files for different data sets of the 
beam-dump data sample. The magnet currents are in Amperes. SMO was off 
during the beam-dump run. 



data sets 


SMI 2 current 


SM3 current 


trigger matrix 


12 


-2800 


-4230 


trigmat.psidump, trigmat.psi2800 


13 


-2040 


-4230 


trigmat.psi2000 


14 


+2040 


+4230 


trigmat.psi2000 


15 


+2800 


+4230 


trigmat.psi2800 



chapter. A total of 82 magnetic tapes of raw data were used, with an average 
of 1.8 Gb of raw data written onto one tape. About 400M events in total were 
recoded during the beam-dump run. 

One data set was distinguished from another by changing magnet-setting 
configurations or trigger-matrix configurations. There were four data sets in the 



beam-dump data sample. The specifications are given in Table |37l . 

The data were taken under two different SM12 settings and two polarities 
in order to reduce possible systematic errors. It was known that the incident 
beam was not perfectly lined-up with the Z-axis of the spectrometer, so the data 
show an up-down asymmetry in the event distribution with respect to the Y = 
plane. The flipping of the magnet polarity thus provided crucial information on 
the measurement of this asymmetry. The changing of SM12 current also changed 
the acceptance of the spectrometer. The consistency on the results obtained 
from the two magnet settings would provide a test, since the physics should be 
acceptance-independent. Each of the four data sets contained about equal amount 
of data. 
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Data set 12 was further divided into data set 12a, 12b, and 12c, and data 
set 15 was divided into 15a and 15b, according to the incident beam angle and 
trigger matrix file. However the beam angle was only determined after the data 
were analyzed, so not till the later chapters will such division be used. 

3.2 Track Reconstruction 

In this section the methodology of track reconstruction applied in the E866 anal- 
ysis is described. For each event, the procedure can be considered as two main 
steps. The first step involved track finding, which was based on the drift-chamber 
hits and muon identification from the proportional-tube signals, and track fitting, 
in which the possible candidates of track segments between Stations 2 and 3 were 
found. In the second step a trace-back procedure was applied to the track candi- 
dates, so that the complete trajectories through the SM12 and SMO to the target 
position were reconstructed. 

3.2.1 Identifying Drift Chamber Hits 

Each drift chamber station consisted of six planes: Y, U, V, and their associ- 
ated prime planes Y',U',V'. When a charged particle traveled through the drift- 
chamber array in one station, correlated signals from different planes were pro- 
duced. The subroutine DCTRIPS searched station 3 and station 2 for the cor- 
related hit patterns. Only the patterns that consisted of at least 4 crossed hits 
whose cross-intersections were very close to a space point were registered. A 
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triplet pattern was defined as having hits in all three views, while a doublet pat- 
tern was defined as having hits in two views only. An associated hit was defined 
as a particle that hits both the prime and unprimed planes. 

3.2.2 Fitting the Tracks 

Once all the valid hit patterns were registered, the next step was to link the reg- 
istered hits from Station 3 and Station 2. The subroutine DCTRAX looped over 
the triplets and doublets in Stations 2 and 3 to construct the track candidates, 
called DC track segments. Several constrains were imposed on the track candi- 
dates: 1) if a doublet in one station was found, it was only allowed to connect to 
a triplet from another station; 2) at least 3 associated hits from Stations 2 and 3 
were required to construct the track segment; 3) the segment was extrapolated to 
Station 4 and was required to fire at least 3 out of 5 planes in the desired location; 
4) the segment was approximately pointing to the target location. For this study, 
very loose cuts were made to confine the segment vectors, and muon identification 
was done in DCTRAX. 

The next step was to link the track segments with the identified hits in 
Station 1. The subroutine WCTRAX required each of the track segments to be 
lined up with a valid hit of Station 1 in the X-Z view (non-bend plane) within 
a vertical band. Only hits within this Station- 1 window were further considered. 
A single bend-plane approximation was used to account for the SM3 momentum 
kick. Once a valid hit was identified in station 1, the entire track was refit into 
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two straight-line segments joint at the bend plane. The fitting routine, FITTIME, 
using all 18 planes of drift chambers to fit the track, and routine SM3 calculated 
the momentum kick and the Z-coordinate of the bend plane. The result of this 
final fit gave the coordinate of the track at the SM3 bend plane, the Y-slopes 
before and after the bend plane, and the X-slope at the intersection point. With 
the knowledge of the SM3 field map, together with the slope information, the 
track momentum at Station 1 was determined. 

3.2.3 Tracing Back through SM12 
3.2.3.1 Energy-Loss Correction 

Prom the SMS bend plane to the target position, the track was reconstructed 
in a routine called PBSWIM. Given the field map of SM12, the coordinates and 
the momentum of the track were reconstructed in the field-map grid step by the 
routine TRACER. During the procedure of tracing back, the effect of energy loss 
in the absorber wall and in the beam dump material was taken into account. The 
lost energy, calculated by an empirical formula 

Eioss = a + 6 * \og{Pin) + c * log(Pi„)' (3.1) 

for each layer of the absorber materials, was added back to the track after TRACER 
had traced through that layer. The coefficients in the formula were determined 
from dedicated Monte Carlo studies. The total energy loss in the beam dump was 
estimated in the same way. For each step inside the dump, a fraction of the total 
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estimated energy loss, proportional to the step size, was added back to the track. 
3.2.3.2 Multiple Scattering 

Due to multiple scattering, it was impossible to trace back to the exact event- 
producing location. So it was assumed that all the events came from a point 
located at one interaction length into the dump. In this case Tjiarget was set to 
86 inches in the E866 coordinate frame. To correct for the effects of multiple 
scattering, a scattering bend-plane approximation was used. After the initial 
traceback, the intercepts of the track at Tjtarget were compared with the beam 
centroid. 



dX -X tar get -^centroidi ('^•^) 

dY ^target ^centroid- (^■^) 

Based on these differences, an angular correction to the track direction at the 
scattering bend plane (located at 'Lscatter) was calculated: 



dOx — dX/[Z tar get ~ Zi scatter) : (3-4) 
dOy = dY / [Z tar get ~ Z scatter)- (3-5) 

After the angular correction was applied, the track was traced again to Z target 
starting from Z scatter- The iteration procedure was repeated until the intercept 
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errors became negligible. The value of Zgcatter was determined by optimizing the 
angular resolution at the target point. 
3.2.3.3 Additional Angle Corrections 

Further Monte Carlo study had shown that the single scattering bend-plane ap- 
proximation actually over-calculated the reconstructed angle, as shown in Figure 

In order to reconstruct the opening angle correctly, an empirical formula, 
the angle deviation expressed as a polynomial function of the reconstructed angle, 
was used to adjust the angle that came out of the initial scattering-plane approx- 
imation. The formula was purely empirical and relied completely on Monte Carlo 
studies, so it was important to test whether these corrections gave back the thrown 
angular distribution for Monte Carlo events after the events were analyzed, even 
though this self-consistency check is only necessary but not sufficient. Figure |3.2| 
shows the reconstructed cos 6 distributions for both magnet settings. Those plots 
were obtained from the thrown-reconstructed events divided by the unsmeared 
acceptances. The function pi x (1 + p2 * cos^ 9) was used to fit the plots to test 
whether there is any systematic bias while reconstructing the cos 6 distributions. 
As a result, a nearly fiat distribution of cos 9 (A = —0.02 ± 0.017) is recovered for 
the 2040Amp data set, and A = 0.02 ± 0.018 for the 2800Amp data set. 

It was also important to test with Drell-Yan data, where we believe we know 
the angular distribution, to search for additional systematic problems not revealed 
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Figure 3.1: The single scattering-plane approximation fails at large angles. The 

reconstructed angles are greater than the thrown angles by about 10% at large 
angles. The figures shown here are taken from the up-going muon tracks. 
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Figure 3.2: Reconstructed cos 6* distributions for both magnet settings. The re- 
constructed distributions recover the thrown (flat) distributions. 
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from Monte Carlo studies alone. The angular distribution study of Drell-Yan 

events as a confidence check will be presented in Chapter 5. 



55 



CHAPTER 4 



DATA ANALYSIS 



Prom the raw-data tapes to the final physics results, several types of computing 
hardware and software were used. The whole procedure can be divided into three 
passes in the analysis, which will be described in this chapter. 

4.1 Pass-1: Fermilab IBM Farm 

The first pass of data analysis was performed on the Fermilab parallel computing 
farms. The computing farm system is a cluster of IBM workstations which can 
distribute the raw data into all the computing nodes and analyze the data simul- 
taneously. The events passing the first-pass analysis were then written on to Data 
Summary Tapes (DSTs) for the second-pass analysis. 

The tasks of the first-pass analysis were mainly to find dimuon events origi- 
nating in the beam dump within the desired mass range. An 18-inch grid SM12 
magnet map was used to trace the tracks in this pass, and very loose aperture and 
target cuts were applied. However it required the events to have two muon tracks 
and the mass of the muon pair has to be greater than 2 GeV. As a consequence 
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only about 5% of the events passed the cuts and were written onto DSTs. 
4.2 Pass-2: Hewlett-Packard Workstation 

The second pass of the data analysis was performed on the Hewlett-Packard Work- 
station located in New Mexico State University. The inputs of this phase were 
DSTs, and the outputs were the ntuple files. The main task of this phase of analy- 
sis was to reconstruct the kinematics of the DST events as accurately as possible. 
A 2-inch grid SM12 magnet map as well as a Y-field map were used to reconstruct 
the events instead of the 18-inch map. No other tighter cuts were applied, but 
many fine-tuning tasks were done in this pass of analysis. 

4.2.1 Determining Tweeks 

The "tweek" is an overall correction factor for the magnetic field strength pro- 
vided by the field map. The field maps provided by the ANL (Argonne National 
Lab) group assumed that the magnets were operated at the preset currents, SM3 
at 4260 Amperes and SM12 at 2800 Amperes, for example. But in reahty the op- 
erating currents were not precisely equal to the desired currents, and there were 
uncertainties in the mapping, therefore it was necessary to apply the corrections 
to the magnet maps for analyzing the events or to generate Monte Carlo events. 

Since the actual currents were not known, this whole subject relied on care- 
ful Monte Carlo studies. There were two unknown quantities to be determined: 
the tweek of SM12 and the tweek of SM3. The two conditions used to deter- 
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Table 4.1: The tweek values of SM12 and SMS for all the data sets. 



Data Set 


SM12(Amp) 


SMS (Amp) 


TWEEK12 


TWEEKS 


12 


-2800 


-4260 


1.006 


0.986 


13 


-2040 


-4260 


1.019 


1.002 


14 


+2040 


+4260 


1.019 


1.002 


15 


+2800 


+4260 


1.006 


0.986 



mine these two quantities were 1) reconstructed J/ip mass and 2) the uniterated 
Z- vertex (ZUNIN). By adjusting the magnet map in the MC event-generating 
phase, it was required that the reconstructed experimental data have the same 
mass and ZUNIN location as the MC reconstructed events; the tweeks used in 
the data event reconstruction were the same as in the MC event generation and 



reconstruction. Figures ^]1| and ^]2| show the ZUNIN and mass peaks of the real 
events and MC events from both magnet settings. The shape of the ZUNIN peaks 
is not symmetric around its central value because of energy loss and multiple scat- 
tering of the muon tracks. The peaks were fitted to a Gaussian using asymmetric 
boundaries, —10 inches to SS inches, in order to locate the peak centroids without 
being affected by the non-Gaussian tails. The mass peaks of the data distributions 
were fitted to a second-order polynomial plus a Gaussian function, while the mass 
distributions from Monte Carlo were fitted to Gaussians since the Monte Carlo 



did not include any background events. As one can see from Figures [4.1| and ^ 
the agreement between Monte Carlo and the experimental data is satisfactory. 
The tweek values are given in Table [4.1| . 
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Figure 4.1: ZUNIN peaks of Data and MC from both magnet settings. 
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Figure 4.2: Mass peaks of Data and MC from both magnet settings. 
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4.2.2 Determining Beam Positions 

During the course of analyzing the raw data, it was found that the beam centroid 
was not steadily fixed at a single location over the entire period of running. The 
moving range of the centroid is greater than three sigma of the beam profile, 
which is 0.14 inches in X and 0.07 inches in Y at the target position, so this was 
due to the beam-line magnet- current fluctuations. A typical reconstructed beam 



centroid distribution during one run is shown in Figure 

Since the angular distribution was affected by the beam angle (a slight offset 
of the beam angle will generate a linear term of cos 9 in the final A extraction) and 
a single fixed target point in the analysis was not able to account for the angle 
fluctuations due to the beam-line movement, the beam centroids at the target 
were determined in a spill-by-spill basis by fitting the beam profile of the raw 
data at the target position. In the first second-pass analysis, the beam centroid 
of each event in each spill was determined by averaging the X and Y coordinates 
of the two muon tracks at the Z = 86 inches plane and then saved the centroid 
distributions (one for X and one for Y) in a temporary histogram. Then these 
distributions were fitted using Gaussian functions to determine the central values. 
These central values were tagged with the run number and the spill number and 
were saved in a 2-D lookup table. This procedure was applied to every spill of 
the raw data. If there were not enough events (the threshold number of events to 
perform the fitting was set to 50) in a spill to perform the fit, the centroid value 
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Figure 4.3: Beam centroids were moving during the data taking. The vertical 
scale gives the Y-position and the horizontal scale gives the X-position in inches 
in the lab frame. 
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Table 4.2: Beam angles of each data set. 



Data set X angle(lE-3) Y angle(lE-3) 



12a 0.12 -1.09 

12b 0.11 -1.05 

12c 0.11 -0.94 

13 0.00 -0.76 

14 -0.02 0.11 
15a 0.10 0.12 
15b 0.05 0.31 



was taken from the previous spill. If the first spill in a run did not have enough 
statistics, the centroid value was set to some default value depending on the data 
set. Then a second second-pass analysis was performed, using the centroid values 
stored in the lookup table, for each spill of data. The beam angles reconstructed 
in this procedure, on a run-by-run basis, were approximately constant within a 
data set. Some runs within a data set had very different angles and had to be 
treated separately. Data set 12 was broken into 12a, 12b, and 12c, and data set 
15 was broken into 15a and 15b. The reconstructed beam angle for each data set 



is shown in Table 4.2. 



4.2.3 Determining Beam Angle 

The determination of the beam angle for each data set involved two phases. In 
the first step we obtained the initial value of the angle recovered by plotting the 
momentum vector from the ntuple. Note that the beam-centroid-fitting procedure 
described in the previous section had to be applied first. Then the second-pass 
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analysis was done again using the initial angle. It was found that after one it- 
eration the reconstructed angles converged within 0.03 mrad. Thus a second 
iteration was not necessary. Figure [4.4] shows the beam angle reconstructed from 



the experimental data set 15b as an example. 

The second step in determining the beam angle was to look at the production- 
(f) distribution. The pro duct ion- (PPHI) distribution was very sensitive to the 
input beam angle. Monte Carlo events were generated to compare with the PPHI 
plots of the real data. The distribution of the real data is expected to be 
isotropical because both the beam and target were not polarized, and in the 
Monte Carlo the PPHI distribution was thrown isotropically. It was required 
that the ratio of the PPHI plots be fiat, so there was no 0-term contribution 
in extracting the angular distribution. Usually small changes in the beam angle 
had to be added to the Monte Carlo in the generation phase to obtain good 
PPHI agreement. However those changes were small compared to the beam angle 
variation of 0.3 mrad, so it was not necessary to repeat the second-pass analysis. 



The ratio of real-data PPHI over MC data PPHI is given in Figure [4.5| . The plot 
was fitted to pi x sin(p2 + 9) + j93. The amplitude shows that the uncertainty of 
the beam angle values used in Monte Carlo is about 2%. 

4.3 Pass-3: PAW 

The third and final pass of the analysis was performed with the Physics Analysis 
Workstation (PAW) program. This program provided a way to present a graphical 
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Figure 4.4: Angle reconstruction plots for data set 15b as an example. The peaks 
are fitted to a Gaussian function to determine the central values. 
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Figure 4.5: Production-^ distribution: experimental data divided by Monte Carlo. 
The plot is fitted to a sine function. The result shows a 2% uncertainty in the 
incident angle. 
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output of the analysis. The input of this program can either be an existing 
histogram file, which is simply a vector, or an ntuple file, which consists of columns 
of vectors. When filtering the final candidates, some additional cuts were applied 
to the events. Target vertex cuts and trigger-bit cuts were imposed. Events 
originating outside a 2"(X)x2"(Y)x200"(Z) range centered at the nominal target 
position were discared. The trigger-bit information TBRAW_CK of each event 
was checked to assure that it satisfied the physics-trigger requirements (satisfies 
Physics Trigger A 1, 2 or 4, see section 2.5.2). Like-sign pairs were discarded by 
the trigger cut at this point. As a result, about 94% of events survived the vertex 
cut and 91% of events survived both vertex and trigger cuts. 

Those events left were then distributed according to their pt and xp values. 
The xp range in this study was from 0.25 to 1.0 binned in intervals of 0.1; pt 
was binned in intervals of 1 GeV. The last bin of xp combined all the data above 
0.85, and the last bin of pt included the data having pt > 3.0 GeV. The mass 
spectrum, from 2.0 GeV to 7.0 GeV plotted in 50-MeV bins, for each xp and pt 
bin was then fitted to a Gaussian plus some background function. The J/V' peak 
was described by the Gaussian function, but the background shape varied as the 
kinematic range changed. Listed below are the functional forms used to do the 
background fitting in this study: 



f{x) — exp{pl -\- p2 * x) , 
67 



(4.1) 



/(x) = pi +p2 * X +p3 * (4.2) 
fix) = pl/{l + {x/p2Y'). (4.3) 

The uncertainties caused by the background function forms are discussed in 
section 4.6. The counts of J/ip^s in each bin of xp and pt were then determined 
by the formula 



COUNTS = N * {bin width/ V2TTa) * exp{-{x - centroidy /2a') (4.4) 

where N, centroid, and a are free parameters to fit. The value + AA^ returned 
from the PAW fitting program provided the population of J/^'s and the statistical 



uncertainty in that bin. Tables |4.3| to |4.6| give the approximate number of J/ijjs 
in each bin for each data set. 



Figure |4.6| shows the reconstructed spectra of some kinematic variables. The 



mass spectrum, on the upper- left, presents all the dimuon pairs recorded during 
the beam-dump run with masses up to 7.0 GeV. The other three variables, namely 
xp, Pt, and cos 6*, are plotted for the events that satisfied the vertex cuts and the 
trigger cuts and have a mass in the range between 2.5 GeV and 4.0 GeV. The 
purpose of the mass cut was to reduce the contributions from non-J/if) events in 
these variables. One can see from the figure that this data sample contains a large 
collection of J/ip's that extend over a wide kinematic range. 
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Table 4.3: Number of J/'^'s in each bin of data set 12. 



/ pAGeV) 


- 1 


1 - 2 


2 - 3 


> 3 


0.3 


322000 


292000 


78800 


16600 


0.4 


531000 


497000 


137000 


28000 


0.5 


323000 


310000 


84800 


16900 


0.6 


142000 


139000 


36500 


6900 


0.7 


43100 


42000 


10700 


1900 


0.8 


8800 


7800 


2200 


300 


> 0.85 


900 


900 







Table 4.4: Number of J/ijj^s in each bin of data set 13. 



/ PT(GeV) 


- 1 


1 - 2 


2 - 3 


> 3 


0.3 


363000 


336000 


91800 


17800 


0.4 


244000 


232000 


66400 


14000 


0.5 


119000 


116000 


32400 


7100 


0.6 


47500 


45900 


13100 


2600 


0.7 


13400 


13100 


3600 


700 


0.8 


2500 


2800 


700 


100 


> 0.85 


300 


300 
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Table 4.5: Number of J/'^'s in each bin of data set 14. 



/ Pr(GeV) 


- 1 


1 - 2 


2 - 3 


> 3 


0.3 


372000 


355000 


99300 


19100 


0.4 


244000 


243000 


71000 


15200 


0.5 


120000 


120000 


34200 


7300 


0.6 


47600 


47300 


13200 


2800 


0.7 


13300 


13600 


3800 


700 


0.8 


2600 


2700 


800 
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> 0.85 


400 
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Table 4.6: Number of J/ijj^s in each bin of data set 15. 



/ Pr(GeV) 


- 1 


1 - 2 


2 - 3 


> 3 


0.3 


260000 


235000 


62300 


13200 


0.4 


413000 


387000 


107000 


22800 


0.5 


250000 


240000 


65200 


13400 


0.6 


109000 


106000 


27700 


5400 


0.7 


33100 


32100 


8400 


1500 


0.8 


6400 


6300 


1500 


300 


> 0.85 


700 


600 
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Figure 4.6: Reconstructed spectra of some kinematic variables. Mass (2.5-4.0 
GeV), vertex, and trigger cuts have been applied to the xp, Pt and cos^ plots. 
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4.4 Monte Carlo 

To extract J/ if) angular distributions correctly we rely on good knowledge of the 
angular acceptance. Monte Carlo simulation is the standard technique to obtain 
the acceptance. In the Monte Carlo study the experimental apparatus setup was 
programmed as close as possible to the real experiment, and all the apparatus 
input parameters, some of which were physically measurable, were tuned accord- 
ing to the best of our knowledge. However the physics part of the Monte Carlo 
generation, which is the part of real interest, was unknown and relied completely 
on theoretical model calculations. The output of the Monte Carlo was then com- 
pared with the experimental data, assuming that the simulation of the apparatus 
part was reliable and trustworthy. The difference between the Monte Carlo re- 
sults and the experimental measurements was then used to improve the various 
thrown physics kinematic quantities. The whole procedure was an iterative pro- 
cess since the acceptance depended on the thrown distributions, and the thrown 
distributions, usually taken from the experimental data distributions, relied on 
the knowledge of the acceptance. The accepted Monte Carlo events were stored 
in the same format as the raw data, and then were analyzed as the experimental 
data. The final output was stored in the form of ntuple files, hke the real data. 
The acceptance correction of a physical variable to be applied to the experimen- 
tal data was given by the reconstructed MC distribution divided by the thrown 
distribution. 
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In this section we will describe the thrown functions of the physical variables 
and compare them with the experimental data. 

4.4.1 xp Distribution 

xp can be understood as the J/ip longitudinal momentum Pl divided by its max- 
imum kinematically allowed value PL,Max, which is approximately equal to half of 
the square root of the center-of-mass energy S, in the beam-target center-of-mass 
frame. Theoretically the xp differential cross section is of interest because it can 
be calculated based on the knowledge of the parton distributions and some phe- 
nomenological models. In this study we used an empirical formula for the thrown 
xp distribution for J/ip within the range 0.25 < xp < 1.0: 



Xp = — ^ —j= — , (4.5) 

PL,Max vS/2 

da/dxp = P3(l - 0.82xp f-'^ , (4.6) 

P3 = 2.784 - 10.14xir + 17.81a;| - 9.585a;^ . (4.7) 

The third-order polynomial was used to describe better the high-Xi? part. The 
comparison of the Monte Carlo and the real data is shown in Figure |4.7|. At each 



bin of Xp the counts of J/ip^s for the real data were obtained by fitting the mass 
spectrum to a Gaussian plus a exponential background. The agreement is good 
to about 3% for xp < 0.7, which includes 98.2% of all the data. The largest 
discrepancy, however, about 20%, for xp > 0.7 comes from the thrown shape in 
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the Monte Carlo and the background uncertainty in the data. This has no effect 
on our results of da / d cos 6 dxp and da / d cos ddxpdpr since these are differential 
quantities. The effects on da /d cos 9 are also expected to be very small because 
only a few events are at that high xp. 

4.4.2 pt Distribution 

Pt is the transverse momentum of the dimuon pair. The origin of Pt is understood 
as a combination of intrinsic transverse motions of the partons inside the hadrons 
and higher-order QCD processes. Naively pr was expected to be an independent 
variable from xp, which accounts for the longitudinal part of the dimuon pair mo- 
mentum, except at some extreme kinematic ranges where the maximum available 
energy becomes a constraint. However experimentally it was found that < pt > 
was correlated with xp beyond pure acceptance effects. In our Monte Carlo code 
we have the following form for the Pt thrown function: 



da/dpT = Pr/il + {PT/p^{xF)ff, (4.8) 
pO = 1.43 + S.2Sxf - lh2xl + 2.66a;^ + 13.9a;^ - 9.23a;f,. (4.9) 

pO was expressed in a polynomial form of xp. This form attemped to fold in 
the real physical correlations between xp and < pp > s-s well as possible residual 



acceptance effects. Figure ^]8| shows < > as a function oi xp for the experi- 



mental data and Monte Carlo data. The overall integrated pp distribution from 
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Figure 4.7: da/dxp of data over da/dxp of Monte Carlo. The agreement is good 
within 3% lor xp < 0.7, in which contains 98.2% of the data. 
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MC was then compared to the integrated pt distribution of the data; this is shown 



in Figure For each bin of p^, the mass spectrum of the data was fitted to 
a Gaussian plus a second-order polynomial as the background. The number of 
J/ip^s was then calculated from the Gaussian parameters. The Monte Carlo cos^ 
distribution was weighted according to the normalized ratio of Figure |4.9| . 

4.4.3 Angular Thrown Distributions 

There are three independent angles used to specify the dilepton pair production: 
decay-^ (DTHETA), decay-0 (DPHI), and production-0 (PPHI). The DTHETA 
variable was defined as the polar angle in the CoUin-Soper frame, and the DPHI 
variable was defined by the azimuthal angle with DPHI = pointing up with 
respect to the Z axis in the C-S frame. The PPHI gave the azimuthal angle of the 
virtual photon around the lab Z-axis. PPHI = was chosen to be the positive 
X direction of the lab coordinate system and it is a lab-frame variable. All three 
angles in this MC code were thrown as flat distributions. The angular distributions 
obtained from the final ntuple thus gave the angular-acceptance shape directly. 

4.5 Extracting Angular Distributions 

In this section the technique for obtaining the polarization parameter A is de- 
scribed. 
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Figure 4.8: < pr > vs xp- Solid line: Data. Dashed line: Monte Carlo. 
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Figure 4.9: da/dpr of data over da/dpr of Monte Carlo. The plot is fitted to a 
zeroth-order polynomial. The result of the fit is the normalization constant used 
in the cos 9 weighting. 
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4.5.1 General Procedure 

da /d cos 9 dprdxp of J/ip was determined by taking the accepted events divided by 
the acceptance, which was obtained by taking the Monte Carlo cos 9 distribution 
divided by the flat thrown distribution. The statistical uncertainty of the accep- 
tance was small compared to the statistical uncertainty of the data. Before being 
applied to the data, the acceptance curves were weighted by the pr distribution 
of the real data. The uncertainties associated with the px thrown function were 
thus expected to be minimized. In addition, since the cos^ acceptances in this 
study were calculated in small xp and pt bins, the uncertainties due to the shape 
of the thrown xp and pr distributions should be reduced. The acceptance of J/ip 
as a function of cos 6' for each bin of xp and pt of the Monte Carlo simulation is 
given in Appendix A. 

To extract the counts of J/ip, the fitting procedure described in section 4.3 
was applied to each bin of xp, Pt, and cos 9, with bin width being 0.1, 1.0 GeV, and 
0.1, respectively. At large xp and pr bins, the bin width of cos^ was increased to 
0.2 to give better statistics. The fitting of each histogram is shown in Appendix 
B. The number of counts, calculated according to Equation 4.4, in each bin of 
cos^ was plotted to obtain the accepted cos^ distribution. Each accepted cos^ 
distribution was then divided by the corresponding acceptance, and the shape of 
da /d cos 9 dpTdxp of J/%1) was thus obtained. The cos^ distribution for each bin 
of a;i? and pt was then fitted according to equation (1.1): 
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da/dcos6 ^ 1 + Acos^ 6. 

The fitting plots of da /d cos Odprdxp are sfiown in Appendix C. 
4.5.2 Combined Data Set 

As mentioned in section 4.2, there were seven data sets differing by the incident 
beam angle, the trigger matrix, or the magnet settings. To obtain the final results, 
those sets were combined according to the SM12 currents. Data sets 12a, 12b, 
12c, 15a, and 15b together formed the "SM12 = 2800" set and data sets 13 and 14 
formed the "SM12 = 2040" set. The events from different experimental data sets 
were directly added together, and the average cos 9 acceptance was calculated by 

< Acceptance >= x aj(cos6')/A^i(cos^) (4-10) 

i 

with 

(4-11) 

i 

where fi is the fraction of the accepted events of data set i among all accepted 
events, and ai(cos6')/A'j(cos6') is the cos^ acceptance distribution of data set i. 
The combined results were obtained by dividing the sum of the accepted events 
by the average acceptance. The results for the two different magnet settings were 
obtained separately. 

In the next chapter and throughout Appendix A to C, the results are pre- 
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sented for the two magnet settings. The agreement between the two sets of results 
provides an important check that our results are not affected by apparatus effects. 
Combined results from the two magnet settings are derived to compare with re- 
sults from other fixed-target experiments. 

4.6 Uncertainties 

In this section two types of errors are discussed: statistical uncertainties and 
systematic uncertainties. These different sources of errors are summarized in 
tables. 

4.6.1 Statistical Uncertainty 

Statistical errors in determining A are the direct results of statistical uncertainties 
in the counts of J/il^. The statistical errors of the acceptance are supressed by 
outnumbering Monte Carlo events over the experimental data. The statistical 
errors of A were obtained by including only the statistical uncertainties of the J / ip 



counting in the A fitting. The A's and the statistical errors are shown in Table 5.1 



and Table 5.2 



4.6.2 Systematical Errors from Analysis and MC Inputs 

To estimate the errors caused by the uncertainties of the magnetic fields, the tweek 
values used in the analysis were varied by 1% in the analysis and Monte Carlo. 
The 1% uncertainty is a reasonable upper limit for the magnetic-field strength 
because the tweek values were tightly constrained by the reconstructed J/ ip mass 
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and the uniterated Z-vertex position. It was found that the A values changed by 
±0.01 as the field strength of SM12 was adjusted, and varied by ±0.01 as SM3 
was changed. 

The incident beam angle has a strong effect on the cos 9 distribution. Al- 
though the value of the beam angle was tuned to remove any asymmetry in the 
ratio plots of decay 9 and production between experimental data and Monte 
Carlo, the precise incident angle is actually unknown. To study the effect of this 
uncertainty, the beam angle used in the analysis and Monte Carlo was varied by 
±0.0002, which is twice the beam angle-spread sigma. The results showed that 
the A values were changed by ±0.02 for large xp and by ±0.04 for small xp- 

The target position in the X-Y plane used in the Monte Carlo was an average 
value determined from the data. In data reconstruction the beam center was 
calculated for each spill. To study this uncertainty, the beam centroid of Monte 
Carlo events was moved by ±0.1 inch in both the X and Y direction. The circle of 
the 0.1 inch confinement was determined by the data distribution. The net effect 
on the A values is ±0.02. 

The pt dependence of the cos 9 distribution and the impact of the px thrown 
function on the cos 9 acceptance was also studied. The cos^ acceptance was 
calculated by weighting the pt distribution according to the real data. Another 
calculation was performed without weighting to the real data. The < pt > with 
and without weighting differed by 5%, which is compatible with the uncertainties 
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Table 4.7: Summary of the systematic errors from all the sources. 



xf SM12 SMS angle centroid < px > fit limit background overall 



0.3 


0.01 


0.01 


0.04 


0.02 


0.06 


0.05 


0.06 


0.109 


0.4 


0.01 


0.01 


0.04 


0.02 


0.06 


0.04 


0.05 


0.099 


0.5 


0.01 


0.01 


0.02 


0.02 


0.06 


0.04 


0.03 


0.084 


0.6 


0.01 


0.01 


0.02 


0.02 


0.06 


0.03 


0.02 


0.077 


0.7 


0.01 


0.01 


0.02 


0.02 


0.06 


0.03 


0.02 


0.077 


0.8 


0.01 


0.01 


0.02 


0.02 


0.06 


0.03 


0.02 


0.077 


0.9 


0.01 


0.01 


0.02 


0.02 


0.06 


0.03 


0.02 


0.077 



in extracting the pO parameter in the px thrown formula. It was then found that 
within this variation, the A values moved by ±0.06. 

4.6.3 Systematical Errors from Peak Fitting 

Other contributions to the systematic errors come from the J/ ip peak fitting pro- 
cess. The A values were found to change slightly with different choices of the 
fitting limits. With the same background function, changing the fitting limit was 
equivalent to changing the continuum shape within the uncertainties. It was found 
that this contributed a ±0.03 uncertainty to the systematic errors. 

The selection of the background function form also produced systematic er- 
rors. Three different functions, equations 4.1, 4.2, 4.3, were used to fit the con- 
tinuum distribution. The uncertainty from different background functions on A is 
±0.02 for large xp and is about ±0.05 for small xp- 



Table |4.7| gives a summary of the systematic errors from all the contributions. 
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CHAPTER 5 



RESULTS 



The angular distributions of the J/V' decay in the dimuon channel have been 
measured for the process p + Cu — > J/tl) -\- X using an 800 GeV proton beam. In 
this chapter the results are presented, along with comparison to the theoretical 
predictions and the results from other experiments. 

5.1 Drell-Yan Angular Distribution 

In the same data sample used in this study, about 200K dimuon pairs with mass 
ranging from 4.0 GeV to 7.5 GeV were also recovered, of which most are Drell- 
Yan events. This data sample is of interest because the target was copper and 
no angular-distribution measurements had ever been published for the proton- 
induced Drell-Yan process in this mass range. Though the amount of data sample 
was not enough to study the polarization as a function of xp, it was still useful 
to examine the overall polarization parameter A, which was expected to be equal 
to unity based on the standard Drell-Yan production mechanism, as a cross check 
for the J/%1) angular distribution results. Below the procedures are described in 
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more detail, and the results are presented 
5.1.1 Random Background 

Random muon pairs were the most significant background contamination in the 
mass range of interest. The definition of a random pair is that two opposite-sign 
muons, which were produced independently by pion decay or other processes, 
coincidentally fired the trigger system and appeared to be a valid target dimuon 
event. The random pair distribution could not be measured directly, because the 
pairs were indistinguishable from the real Drell-Yan dimuons in the spectrometer; 
they were however simulated from the like-sign event distribution by changing 
the sign of the Y-momentum of one of the like-sign tracks to calculate other 
kinematic variables of the pair, based on the assumption that the probability to 
form a random pair is the same as to form a like-sign pair given the first muon 
track. Single-muon trigger rates were used to normalize the ratio of like-sign pairs 
and the randoms; the like-sign pairs were expected to have the same single-muon 



trigger rate as the random pair would have. Figure |5.1| shows some kinematic- 
variable distributions of the random pairs. 

5.1.2 Random Subtraction and Results 



Prom Figure ^.2| it was understood that the random subtraction was important in 
order to extract the Drell-Yan angular distribution, even for pairs of mass greater 
than 5.0 GeV. The normalization factor for the random events was determined 
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Figure 5.1: Some kinematic variables of the random pairs. The mass range of the 
xp, Pt: and cos^ plots is above 2 GeV. 
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Figure 5.2: 0086* of the random pairs. Solid line: pair mass > 4 GeV. Dashed line: 
pair mass > 4.5 GeV. Dotted line: pair mass > 5 GeV. 
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by matching the number of hke-sign pairs from the random ntuple file and from 
the data, since the random events were generated according to the amount of 
hke-sign pairs of the data. Then the angular acceptance of the Drell-Yan events 
was obtained by a dedicated Monte Carlo run. The cos 9 distribution of the Drell- 
Yan data, after subtracting out the randoms, was corrected for acceptance and the 
angular distribution was obtained. The angular distribution of Drell-Yan events in 



the mass range of 4.0 GeV to 7.0 GeV is shown in Figure |5.3| . The rise at the edges 
is understood as a resolution problem and the same effect has been reproduced by 
broadening the cos 9 resolution. Also, the random pairs show a very strong rise at 
large cos 6*; a slight mismatch in the normalization can result in the same effect. 
For these reasons, for the Drell-Yan data the angular distribution was fitted in 
the range of —0.7 < cos9 < 0.7 in which the systematic uncertainties are best 
handled. A result of A = 0.98 ± 0.04 was obtained. This is consistent with 1.0 
as predicted. This provided a confidence check for the J/ip angular distribution 
presented in next section. 

5.2 J/ifj Angular Distribution Results 

In this section the J/ip angular distribution results are presented. The mea- 
surements were performed under two different magnet configurations, therefore 
different acceptances, in order to minimize the systematic bias. They are effec- 
tively two independent measurements. The polarization parameter A in bins of pt 
and xp is presented in Tables ^TT] and |^, and plotted in Figure |^ . In the figures 
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Figure 5.3: cos^^ of the Drell-Yan pairs with mass ranging between 4 GeV and 7 
GeV. Solid line: After random subtraction. Dashed line: Before random subtrac- 
tion. A A value of 0.98 ± 0.04 is obtained after correcting for the random pairs. 
The rise at the edges is due to resolution effects. 
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Tabic 5.1: A ill Xp and pp bins with statistical errors only- This tabic is for the 

SM12 = 2010 (lata. 

Xp <pt < I 1 <pt <2 2 <pt <3 3 < pr 

0.25 - 0.35 0.153 ±0.037 0.057 ±0.024 0.093 ± 0.026 0.124 ±0.049 

0.35 - 0.45 0.218 ±0.031 0.015 ± 0.019 0.095 ± 0.026 0.141 ±0.056 

0.45 - 0.55 0.146 ±0.023 0.035 ± 0.017 0.101 ± 0.025 -0.052 ± 0.049 

0.55 - 0.65 0.151 ±0.039 -0.013 ± 0.027 0.072 ± 0.041 -0.06 ± 0.08 

0.65 - 0.75 0.111 ±0.070 -0.211 ± 0.046 0.023 ± 0.093 -0.43 ±0.13 

0.75 - 0.85 -0.17 ±0.15 -0.22 ± 0.09 -0.14 ±0.30 
> 0.85 -0.44 ±0.42 



Table 5.2: A in a;^ and pt bins with statistical errors only. This table is for the 
SM12 = 2800 data. 



Xp 


<Pt <l 


1 < Pt < 2 


2<pt <S 


3 < Pt 


0.25 - 0.35 


0.189 ±0.044 


0.135 ±0.023 


0.120 ±0.034 


0.060 ±0.064 


0.35 - 0.45 


0.162 ±0.029 


0.095 ±0.015 


0.170 ±0.021 


0.170 ±0.045 


0.45 - 0.55 


0.115 ±0.022 


0.051 ±0.013 


0.153 ±0.020 


0.057 ±0.040 


0.55 - 0.65 


0.018 ±0.028 


-0.053 ± 0.020 


-0.026 ±0.028 


0.01 ± 0.06 


0.65 - 0.75 


-0.032 ± 0.049 


-0.174 ±0.033 


-0.167 ±0.059 


-0.09 ±0.12 


0.75 - 0.85 


-0.25 ±0.10 


-0.09 ±0.10 


-0.21 ±0.19 




> 0.85 


-0.51 ±0.54 









only the statistical errors are shown. The results from the two measurements are 
in agreement with each other. This provides a confirmation that the results are 
not affected by the specific apparatus settings. 

5.3 Comparison 

In order to compare with other experiments, A values inxp and pt bins from E866 
measurements were combined to obtain X in xp bins and the overall integrated A 
using the following the relations: 
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Figure 5.4: J/ijj polarization parameter X in xp and pt bins. The errors shown 
here are statistical only. 
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Table 5.3: A in xp bins with statistical errors only. 



Xp 


A(xf)(SM12 = 2040) 


A(xf)(SM12 = 2800) 


Combined 


0.25 - 0.35 


0.092 ±0.015 


0.134 ±0.017 


0.110 ±0.011 


0.35 - 0.45 


0.081 ±0.013 


0.129 ±0.011 


0.109 ±0.008 


0.45 - 0.55 


0.073 ±0.012 


0.086 ±0.009 


0.081 ±0.007 


0.55 - 0.65 


0.041 ±0.019 


-0.026 ±0.014 


-0.002 ±0.011 


0.65 - 0.75 


-0.116 ±0.034 


-0.134 ±0.024 


-0.128 ±0.020 


0.75 - 0.85 


-0.200 ± 0.073 


-0.174 ±0.066 


-0.186 ±0.049 


> 0.85 


-0.44 ±0.42 


-0.51 ±0.54 


-0.47 ±0.33 


All data 


0.065 ±0.007 


0.070 ±0.005 


0.069 ±0.004 



yl = ^^ (5.1) 
Wi = Xja] (5.2) 

where cxj is the statistical error of some measurement Aj, and A is the average 



value of Ai. The results are shown in Table |5.3| . One can see from the table 
that the J/-?/; starts slightly transversely polarized at small Xp^ then eventually 
becomes partially longitudinally polarized as increases toward unity. The A 



values from both magnet settings versus xp are plotted in Figure 



The E866 results are compared with the results of CIP data |[Bii 87[|. Figure 



^ shows the E866 results and data published by the CIP group. Recall that the 
CIP experiment was fixed-target vrN collisions. At large xp^ both experiments 
observe longitudinal polarizations. At smaller xp^ E866 sees small transverse po- 
larization while the CIP group saw no polarization. Since the dominant Feynman 

diagrams are different for pN (mainly g-g fusion) and ttN (it has significant qq 
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Figure 5.5: J/ip polarization parameter A in xp bins. The errors are statistical 
only. 
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Table 5.4: Overall A values from other fixed-target experiments and E866. 



Experiment reaction \^ xp range A 

E537 p + W 15.3 GeV xf>^ -0.115 ±0.061 

E537 TT- + W 15.3 GeV > 0.028 ± 0.004 

E672/706 TT- + Be 31.5 GeV 0.1 < < 0.8 -0.01 ± 0.08 

E771 p + Si 38.8 GeV -0.05 < < 0.25 -0.09 ± 0.12 

E866 p + Cu 38.8 GeV 0.25 < x^. < 1.0 0.069 ± 0.004 ± syst. 



contributions) at small xp-, the differences in the polarization are not unexpected. 
However if the 0.1 systematic errors are included, the E866 results at x^? < 0.5 
are marginally in agreement with no polarization. 

The integrated polarization parameter A obtained by E866 and other previous 
experiments are presented in Table |5.4| for comparison. Recall that E866 uses the 



Collin-Soper frame and the other experiments have used the Gottfried- Jackson 
frame as their reference frame. However the pt in fixed-target experiments is low 
enough that the direct comparison is still sensible. E866 gives A = 0.069 ± 0.004 
integrated over all available data. If the systematic error is included, the E866 
result shows no polarization. The transverse polarization at small Xi? is partially 
cancelled by the longitudinal polarization at large xp- The overall result is in 
agreement with other experiments, and in contrast to the non-relativistic QCD 
calculation. 

The A's were also integrated in 1-GeV px bins for xp < 0.45 and xp > 0.45 
to study the pp dependence, as shown in Figure pT7[ No pp dependence was 
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Figure 5.6: J ftp polarization parameter A in xp bins from the FNAL E866 and 

from Chicago-Iowa-Princeton coUaboration. The error bars on E866 data are 
statistical only; the systematic error is shown in the shadowed band below. 
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identified for either xp region. This suggests that nuclear effects are probably 
not responsible for the polarization observed, since one important cause of the 
broadening in the pr distribution in nuclear targets is the multiple scattering of 
the incoming and outgoing partons with the nuclear media. If some of the nuclear 
effects, e.g. energy loss inside nucleus, are important, one would expect to see 
significant px dependence on the polarization. 
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Figure 5.7: J/ip polarization parameter A in 1-GeV px bins for xp < 0.45 and 
xp > 0.45. The errors are statistical only. 
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CHAPTER 6 



CONCLUSIONS AND FUTURE PROSPECTS 



The angular distribution of J/ip decays in the channel produced in 800 

GeV proton-copper collisions has been measured for xp > 0.25. The polarization 
parameter A is extracted in pt and xp bins for two magnet configurations with 
different acceptances. The data indicate that the J/ip^s are produced with a 
slight transverse polarization at xp < 0.6, which turns to longitudinal a.t xp > 
0.6. This suggests that gluon-gluon fusion, which dominates at small xp, and 
quark-antiquark annihilation, which dominates at large xp, leave J/ip's in different 
polarization states. Another fixed-target experiment | Bii 87 |, using pion beams. 



also showed longitudinal polarization at xp 1. However at smaller xp the 
uncertainties are large and no evidence of polarization is seen in Biino's paper. 
The difference of the results from E866 will provide interesting information on how 
the production mechanism affects the polarization, because in ttN interactions the 
production is dominated by quark-antiquark annihilation while in the pN case the 
production is dominated by gluon-gluon fusion in the range of xp < 0.6. 

It should be mentioned that the J/tp samples collected in this study do not 
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Table 6.1: Ratio of primary proton and secondary pion induced- J/^/; [[Mue 99 



xf ratio(%) 

02 15.1 

0.3 9.3 

0.4 6.2 

0.5 4.4 

0.6 3.4 

0.7 3.0 

0.8 2.9 



purely come from pN interactions. A significant amount of pions were generated 
at the dump by hadronic interactions, and those pions can further interact with 



the beam dump to produce J/^/^'s. A calculation ||Mue 99|| shows that about 10% 



of the J /ip^s come from pion interactions at small xp- This has to be taken into 



account when comparing with the theoretical calculations. Table |6.1| shows the 
estimate of the ratios of the J/ip^s produced by the secondary pions to those 
produced by the primary proton beam in various xp ranges. 

It is also important to keep in mind that a substantial fraction of J/ip^s come 
from decays of the Xc states and ip' decays in addition to direct J/ip's. All the 
processes contribute different amounts of polarization to J/ip. Thus one needs to 
know the relative production cross sections of the various charmonium states to 
interpret the results properly. So far, only production ratios for pion-produced 
charmonium states are available. It is also necessary to know the polarization of 
J/ip's from each process to extract the polarization of direct J/ip decays. A theo- 
retical calculation has been done for vrN collisions using the Color-Singlet Model 
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Van 95 1, but the results do not agree with the pion data. It would be interesting 
to see the predictions of similar calculations for pN collisions. A measurement of 
polarization would be very interesting since the Xc state contribution is absent. 
In this experiment the mass resolution was sacrificed to gain the yield rate and 
angular coverage; the data sample contains only about 1 % of ip'^s and they are 
not resolved from the J/ peak. 

It is interesting to notice that if the J/ip^s are integrated over the entire xp 
range, the transverse polarization at small xp partially cancels the longitudinal 
polarization at large xp, and the overall effect appears to be no polarization if 
the systematic uncertainty is included. Unpolarized J/ip^s were also observed in 
other fixed-target experiments, using either proton or pion beams. 

Nuclear effects may also affect the J/ ip polarization, since the J/ if) may collide 
with other nucleons before it can escape the nucleus. The original polarization 
may thus be supressed or smeared out. To eliminate such an effect, a hydrogen 
target is preferable, at the price of smaller production rate however. It would also 
be interesting to study the nuclear dependence of A to understand the nuclear 
effects on the polarization patterns. 

The large-XiT' behavior is of interest and yet remains mysterious. The po- 
larization is changed to longitudinal. Similar behavior was observed also in the 



pion data and a possible explanation is higher-twist effects ||Van 95||. It is not 



clear however how this mechanism applies to proton-induced data. In this study 
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the xp coverage is only up to 0.95. The statistics are too poor to produce sen- 
sible results for xp > 0.95. Even for the 0.85 < xp < 0.95 bin it is desirable 
to reduce the statistical uncertainty. It would be interesting to know whether A 
actually drops to —1 when xp approaches 1.0 in pN interactions. This might give 
us better understandings of the higher-twist effects. 
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APPENDIX A 

cos^ Acceptance in xp and pt Bins 



The cos 6 acceptance inxp and pt bins for both magnet settings is presented in this 
appendix. In Figures Al to A6 the cos^ acceptance for the "SM12=2040" data 
is plotted, and from Figure A7 to A12 the cos 6* acceptance for the "SM12=2800" 
data is plotted. 
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Figure A.l: cos^ acceptance in xp and pt bins of the 2040Amp data. 
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Figure A. 4: cos^ acceptance in xp and pt bins of the 2040Amp data. 
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Figure A. 5: cos^ acceptance in xp and pt bins of the 2040Amp data. 
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Figure A. 6: cos^ acceptance in xp and pt bins of the 2040Amp data. 
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Figure A. 7: cos^ acceptance in xp and pt bins of the 2800Amp data. 
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Figure A. 8: cos^ acceptance in xp and pt bins of the 2800Amp data. 
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Figure A. 9: cos^ acceptance in xp and pt bins of the 2800Amp data. 
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Figure A. 10: cos^ acceptance in xp and pt bins of the 2800Amp data. 
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Figure A. 11: cos^ acceptance in xp and pt bins of the 2800Amp data. 
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Figure A. 12: cos^ acceptance in xp and pt bins of the 2800Amp data. 
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APPENDIX B 
Fitting the J/ijj Peaks 



In this appendix the fitting of the J/%1) peaks is presented. The di-muon mass 
spectrum from each bin of cos^, pr, and xp, are fitted to a Gaussian plus a 
background function with J/^i/^'s fitted to a Gaussian shape. The count of J/ip^s 
in each bin was then calculated according to the output parameters of the Gaussian 
fit. Figure Bl to B42 show the fittings from the "SMI 2 = 2040 Amp" data set, 
and Figure B43 to B84 show the fittings from the "SM12 = 2800 Amp" data set. 

These plots have been removed from this version of the thesis. A complete 
copy of this Appendix is available from the author, or by contacting Department 
of Physics, New Mexico State University, Box 3D, Las Cruces, NM 88003. 
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APPENDIX C 
Fitting the A's 



After the counts of the J/^l) peaks for each bin of cos^, xp-, and pr, were deter- 
mined, the number of counts was plotted versus cos^. This gave us the accepted 
cos 6* distributions for J/ip in bins of xp and pt- Then those accepted cos^^ dis- 
tributions were divided by the acceptance curves shown in Appendix A to ob- 
tain the true distributions. The corrected cos 6* distributions were then fitted to 
1 -|- Acos^^ times an arbitrary normahzation constant. Those distributions and 
the fits are presented in this appendix. In Figures CI to C6 are the curves for 
the "SM12=2040" data set, and from Figures C7 to C12 are the curves for the 
"SM12=2800" data set. 
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Figure C.l: The corrected cos^ distributions and the polarization parameter A. 
Upper left: < pt < 1 and 0.25 < xp < 0.35. Upper right: < pt < 1 and 
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Figure C.2: The corrected cos^ distributions and the polarization parameter A. 
Upper left: < < 1 and 0.65 < xp < 0.75. Upper right: < pt < 1 and 
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Figure C.3: The corrected cos^^ distributions and the polarization parameter A. 
Upper left: 1 < pt < 2 and 0.35 < xp < 0.45. Upper right: I < pr < 2 and 
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Figure C.4: The corrected cos^ distributions and the polarization parameter A. 
Upper left: 1 < pt < 2 and 0.75 < xp < 0.85. Upper right: 2 < < 3 and 
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Figure C.5: The corrected cos^ distributions and the polarization parameter A. 
Upper left: 2 < < 3 and 0.55 < xp < 0.65. Upper right: 2 < < 3 and 
0.65 < xf < 0.75. Lower left: 2 < pr < 3 and 0.75 < xp < 0.85. Lower right: 
3 < pt and 0.25 < xp < 0.35. SM12=2040. 
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Figure C.6: The corrected cos^^ distributions and the polarization parameter A. 
Upper left: 3 < px and 0.35 < xp < 0.45. Upper right: 3 < and 0.45 < 
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Figure C.7: The corrected cos^ distributions and the polarization parameter A. 
Upper left: < Pt < 1 and 0.25 < Xp < 0.35. Upper right: < Pt < 1 and 
0.35 < xp < 0.45. Lower left: < pr < 1 and 0.45 < xp < 0.55. Lower right: 
< pt < 1 and 0.55 < xp < 0.65. SM12=2800. 
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Figure C.8: The corrected cos^ distributions and the polarization parameter A. 
Upper left: < < 1 and 0.65 < xp < 0.75. Upper right: < pt < 1 and 
0.75 < xp < 0.85. Lower left: < pt < 1 and 0.85 < xp- Lower right: 1 < pr < 2 
and 0.25 < xp < 0.35. SM12=2800. 
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Figure C.9: The corrected cos^^ distributions and the polarization parameter A. 
Upper left: 1 < pt < 2 and 0.35 < xp < 0.45. Upper right: 1 < pt < 2 and 
0.45 < xp < 0.55. Lower left: 1 < pr < 2 and 0.55 < xp < 0.65. Lower right: 
1< Pt < 2 and 0.65 < xp < 0.75. SM12=2800. 
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Figure CIO: The corrected cos^^ distributions and the polarization parameter A. 
Upper left: 1 < < 2 and 0.75 < xp < 0.85. Upper right: 2 < < 3 and 
0.25 < xf < 0.35. Lower left: 2 < < 3 and 0.35 < xp < 0.45. Lower right: 
2 < pt < 3 and 0.45 < xp < 0.55. SM12=2800. 
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Figure C.ll: The corrected cos^ distributions and the polarization parameter A. 
Upper left: 2 < pt < 3 and 0.55 < xp < 0.65. Upper right: 2 < < 3 and 
0.65 < xf < 0.75. Lower left: 2 < pr < S and 0.75 < xp < 0.85. Lower right: 
3<Pt and 0.25 < xp < 0.35. SM12=2800. 
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Figure C.12: The corrected cos^ distributions and the polarization parameter A. 
Upper left: 3 < and 0.35 < xp < 0.45. Upper right: 3 < px and 0.45 < 
xp < 0.55. Lower left: 3 < pt and 0.55 < xp < 0.65. Lower right: 3 < pt and 
0.65 <xf< 0.75. SM12=2800. 
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